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ABSTRACT 

A  method  of  designing  airfoils  in  cascade  by  means  of 
conformal  transformations  is  discussed.  Design  of  these  air- 
foils is  regulated  by  five  independent  input  parameters,  with 
solutions  obtained  by  digital  computer.,   A  large  number  of 
cascades  generated  by  this  method  are  compared.   To  evaluate 
the  limits  of  performance,  a  parameter  to  indicate  the  ten- 
dency toward  flow  separation  is  introduced,  with  a  limiting 
value  established  and  verified.  Proper  solidity  is  shown  to 
be  of  great  importance  in  achieving  low  values  of  this  separ- 
ation parameter.  The  value  of  proper  solidity  for  a  given 
blade  thickness  is  shown  to  be  relatively  independent  of 
fairing  shape.  To  increase  performance,  reducing  blade  thick- 
ness with  a  corresponding  increase  in  solidity  is  shown  to 
yield  much  greater  improvement  than  changes  in  fairing  shape. 
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1.0  Introduction 

This  study  is  an  extension  of  a  line  of  investigation 
originally  conducted  by  Professor  Theodore  H.  Gawain  of  the 
U.  S.  Naval  Postgraduate  School.  The  original  investigation, 
sponsored  by  the  Aerojet  General  Corporation,  was  specifically 
aimed  at  increasing  the  pressure  rise  per  stage  of  axial  flow 
pumps.  The  intended  application  was  to  reduce  the  size  of 
propellant  pumps  for  large  liquid  fuel  rockets.  The  require- 
ments were  for  maximum  pressure  rise  per  stage,  regardless 
of  efficiency.  This,  in  turn,  demanded  positive  and  accurate 
control  of  the  pressure  distribution  around  each  blade,  with 
particular  emphasis  on  the  adverse  pressure  gradient  experi- 
enced by  the  boundary  layer. 

Flow  of  a  real  fluid  through  an  axial  flow  compressor  or 
pump  involves  some  formidable  problems.  Besides  the  effects 
of  viscosity  and  turbulence;  tip  losses,  centrifugal  force 
fields,  and  other  effects  complicate  the  problem  considerably. 
As  is  usually  the  case  in  engineering  however,  a  simplifica- 
tion of  the  problem  may  yield  useful  results.  If  the  stream 
surfaces  of  the  flow  are  assumed  to  be  concentric  cylinders, 
then  a  cylindrical  cross  section  of  the  flow  may  be  repre- 
sented by  an  infinite  linear  cascade  with  two-dimensional  flow. 
If  all  velocities  are  considered  relative  to  the  blade  row, 
then  the  solution  for  the  rotor  and  the  stator  have  identical 
mathematical  form.   Following  the  usual  aerodynamic  practice, 
viscosity  may  be  neglected  outside  the  boundary  layer.  This 


leaves  only  turbulence  as  a  factor  preventing  the  use  of  poten- 
tial flow  considerations.   In  order  to  use  the  powerful  analy- 
tic methods  of  potential  flow,  turbulence  has  been  neglected. 
This  does  not  mean  however,  that  the  data  presented  is  invalid. 
Actually,  the  potential  flow  assumption  is  used  to  obtain  the 
pressure  distribution  only.  The  relation  between  pressure 
distribution  and  separation  was  obtained  from  experimental 
data  from  actual  turbomachines,  and  must  therefore  include  the 
effects  of  turbulence. 

From  the  foregoing  considerations,  it  may  be  seen  that 
the  study  of  cascade  design  has  application  to  the  whole  field 
of  axial  flow  turbomachinery.   The  problem  of  cascade  design 
is  much  more  difficult  than  the  problem  of  isolated  airfoil 
design.  In  the  early  days  of  airfoil  design,  much  experimen- 
tal data  was  obtained  on  a  bewildering  variety  of  shapes. 
However,  there  was  no  reliable  method  of  predicting  the  per- 
formance of  an  airfoil  before  it  was  tested.  It  was  only 
after  the  NACA  (now  NASA)  discovered  how  airfoil  performance 
depended  on  shape,  specifically  on  thickness,  camber  and  fair- 
ing distribution,  that  any  real  progress  was  made  in  the  field 
of  airfoil  design. 

The  present  state  of  the  field  of  cascade  design  is  that 
there  is  a  large  volume  of  experimental  data  available  on  a 
bewildering  variety  of  cascades.  Many  efforts  have  been  made 
to  extrapolate  the  theory  of  isolated  airfoils  so  that  it  will 
apply  to  cascades.   However,  at  present  there  is  no  reliable 


theory  which  covers  the  field  of  cascade  performance  in  the 
way  that  NASA's  theory  covers  the  field  of  airfoil  performance. 
The  parameters  of  thickness,  camber,  and  fairing  distribution 
are  still  involved,  plus  a  few  more.  The  relations  and  in- 
teractions between  these  parameters  are  so  complex,  that  it 
is  very  difficult  to  consider  any  aspect  of  cascade  design 
separately. 

When  considering  cascades  as  complete  entities,  the  inter 
action  effects  are  all  present,  and  need  not  be  corrected  for. 
However,  experimental  efforts  in  this  direction  have  been  nec- 
essarily limited  in  scope,  due  to  the  effort  and  expense  in- 
volved in  constructing  and  testing  large  numbers  of  cascades. 
Theoretical  methods,  such  as  the  one  outlined  in  this  investi- 
gation, provide  the  means  of  considering  complete  cascades 
over  a  wide  range  of  design  parameters.  This  study  was  under- 
taken from  the  standpoint  of  discovering  areas  of  possible 
improvement  in  compressor  performance* 

Most  of  the  efforts  of  the  NASA  in  the  field  of  compre- 
sor  design  have  been  for  maximum  efficiency  with  moderate 
output.   In  some  cases,  it  may  be  required  to  attain  maximum 
output,  with  perhaps  moderate  efficiency.  The  gains  to  be 
expected  from  these  considerations  are  not  of  large  order  how- 
ever, as  a  cascade  develops  its  best  efficiency  not  far  from 
the  point  of  stalling  the  blades. 

The  basic  problem,  then,  is  composed  of  two  parts.   First 
it  is  necessary  to  devise  a  method  for  generating  blade  shapes 


with  the  corresponding  pressure  distributionso  For  this  pur- 
pose, the  powerful  methods  of  complex  variables  and  conformal 
transformations  are  utilized.  The  method  developed  by  Pro- 
fessor Gawain  generates  infinite  cascades  of  airfoil  shapes, 
based  on  five  arbitrary  parameters.  Generally,  the  transfor- 
mations are  similar  to  the  Joukowski  transformations  with 
additional  requirements. 

The  second  part  of  the  problem  is  to  establish  a  correla- 
tion between  the  pressure  distribution  around  each  blade  and 
the  occurence  of  separation.  This  problem  must,  of  course, 
depend  heavily  on  experimental  data.  Fortunately,  work  has 
been  done  in  this  area  and  some  practical  working  limits  have 
been  established.  See  Reference  1.  Based  on  this  data,  the 
simple  parameter  (  Vmax/V2)^  has  been  chosen  as  an  indication 
of  the  tendency  toward  separation. 

The  author  wishes  to  express  his  sincere  appreciation  to 
Professor  Theodore  H.  Gawain  of  the  U,  S.  Naval  Postgraduate 
School  for  his  assistance,  guidance,  and  encouragement  during 
the  course  of  this  investigation. 


2.0  The  Method  of  Cascade  Design 

2,1  Introduction 

Methods  of  designing  airfoils  by  purely  mathematical  means 
have  been  known  since  the  early  days  of  aviation,   Joukowski 
was  the  most  notable  pioneer  in  this  field.  He,  and  others, 
noted  the  possibility  of  designing  cascades  of  airfoils  by 
similar  methods.  Since  the  early  applications  of  cascade  de- 
sign were  to  multi-winged  airplanes,  not  much  work  was  done  in 
this  field.  As  pointed  out  earlier,  the  application  to  turbo- 
machinery  makes  this  field  of  considerable  interest  today. 

In  considering  the  flow  through  cascades,  the  shape  and 
position  of  any  individual  blade  will  affect  the  flow  around 
that  blade.  In  addition,  the  shape  and  position  of  all  the 
other  blades  will  also  affect  the  flow  around  that  same  blade 
at  the  same  time.   Therefore,  more  parameters  must  be  used  in 
cascade  design  than  are  used  in  isolated  airfoil  design.  In 
addition  to  blade  shape,  the  relation  of  the  blades  to  each 
other,  and  their  relation  to  the  flow  must  all  be  described. 
The  method  developed  by  Professor  Gawain  uses  five  independent 
parameters  to  describe  an  individual  cascade. 

The  solution  for  blade  shape  and  velocity  distribution 
corresponding  to  chosen  values  of  the  five  parameters  is  found 
by  using  relations  of  potential  flow  and  conformal  mapping  in 
three  separate  planes.  These  have  been  termed  the  near  circle 
plane,  the  circle  plane,  and  the  cascade  plane  respectively. 


In  reading  through  the  following  sections,  which  describe 
the  general  method  of  generating  cascades,  reference  may  be 
made  to  Figure  5>  which  illustrates  a  typical  contour  in  each 
of  the  planes  referred  to  above.  A  summary  of  the  mathematical 
development  of  the  method  is  given  in  Appendix  A.  A  more  de- 
tailed analysis  may  be  found  in  Reference  2. 

The  solution  of  these  equations,  many  of  which  are  trans- 
cendental, is  an  impossible  task  by  hand  calculation.  The 
entire  problem  has  been  arranged  in  FORTRAN  language  for  solu- 
tion by  a  digital  computer.  The  CDC  1604  computer  at  the  U.  S. 
Naval  Postgraduate  School  was  used  for  this  investigation. 
The  basic  computer  program,  along  with  later  modifications  is 
included  in  Appendix  B. 

The  following  sections  describe  the  general  method  of 
generating  cascades  from  five  input  parameters.  Each  parameter 
is  described  in  detail  as  it  enters  into  the  calculations. 
The  parameters  are: 

P     =  Mean  flow  direction 

30    m   Turning  parameter 

O*    ■  Equivalent  flat  plate  solidity 

£*  »=  Thickness  parameter 

/<   =  Shape  or  symmetry  parameter 

2.2  The  Flow  Parameters 

Flow  through  a  cascade  is  usually  described  in  terms  of 
two  velocities.  The  first  velocity  is  that  far  enough  up- 
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stream  of  the  cascade  so  that  local  effects  due  to  the  cascade 
are  not  evident.  This  is  called  the  inlet  velocity,  and  the 
angle  it  makes  with  the  reference  axis  is  called  /3L.     Similarly, 
the  velocity  far  enough  downstream  so  that  local  effects  have 
damped  out  is  called  the  outlet  velocity,  and  its  angle  with 
the  reference  axis  is  called  ($2°     The  magnitudes  of  these  two 
velocities  are  not  independent,  since  the  laws  of  continuity 
require  that  the  component  perpendicular  to  the  cascade  axis 
be  equal  for  the  two  cases.  This  concept  allows  velocities  to 
be  conveniently  nondimensionalized,  so  that  only  angles  are 
important . 

Figure  1  is  an  illustration  of  two  typical  cascades.   An 
appropriate  velocity  vector  diagram  is  shown  for  each.  Note 
that  only  two  parameters  are  necessary  to  describe  the  complete 
diagram.  Instead  of  the  angles  (3±   and  (3z   it  was  more  conven- 
ient to  use  two  different  parameters.  Consider  the  following 
expressions: 

Tan  (3    -  £(Tan  fa   +  Tan  (32) 
B0  =  (Tan  (\   -Tan  /82) 

The  first  expression  defines  the  angle  /?,  or  mean  flow  direc- 
tion, which  is  the  first  input  parameter.  The  second  expres- 
sion defines  the  turning  parameter  BQ,  which  is  the  second 
input  parameter.  The  relationship  between  all  these  quantities 
is  shown  in  the  diagram  on  Figure  1. 

The  starting  point  of  the  calculations  is  the  near  circle 
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plane.  In  this  plane,  the  complex  potential  for  a  series  of 
doublets  spaced  at  intervals  of  2  Tf  along  the  y  axis  is  com- 
bined with  the  potential  for  uniform  flow  at  the  angle  /?  with 
the  x  axis.  The  velocity  of  this  flow  is  set  at  unity,  which 
then  becomes  the  reference  for  all  the  nondimensional  veloci- 
ties derived  later.  If  only  one  doublet  were  used,  the  resul- 
ting potential  would  describe  the  flow  without  circulation 
around  a  single  circular  cylinder.  In  this  case,  with  an  in- 
finite series  of  doublets  along  the  y  axis,  the  resulting 
potential  describes  the  flow  without  circulation  around  an 
infinite  series  of  nearly  circular  bodies  spaced  along  the  y 
axis.  These  bodies  are  not  actually  circles  because  the  pre- 
sence of  the  infinite  stack  of  doublets  distorts  the  flow 
around  any  individual  doublet,  such  that  the  flow  contour  is 
a  slightly  flattened  circle.  The  amount  of  distortion  is  re- 
lated to  the  strength  of  the  doublets  and  the  spacing  between 
them.  In  this  case,  the  spacing  between  doublets  is  fixed  at 
2  Tf   ,  so  doublet  strength  is  what  determines  the  relation 
between  spacing  and  size  of  the  bodies  in  the  near  circle  plane. 

2.3  The  Solidity  Parameter 

The  strength  of  the  doublets  is  fixed  by  the  input  para- 
meter GT  i  or  the  equivalent  flat  plate  solidity.  This  para- 
meter is  closely  related  to,  but  not  identical  with, the  ordinary 
solidity  c/S,  or  the  chord-pitch  ratio  of  the  blade  row.  The 
distinction  can  best  be  illustrated  by  an  example.  Generally, 


t 


any  row  of  arbitrary  blades  of  solidity  c/S  may  be  transformed 
conformally  into  a  related  row  of  flat  plate  vanes.  The  chord- 
pitch  ratio  of  the  latter  is  then  denoted  by  CT  .  It  follows 
that  if  two  arbitrary  cascades  have  the  same  value  of  0"  , 
then  either  may  be  conformally  transformed  into  the  other.  But 
if  the  value  of  (H  is  different  for  each,  then  no  such  trans- 
formation is  possible.   (T"  then,  may  be  regarded  as  a  kind 
of  generalized  solidity,  of  perhaps  deeper  aerodynamic  signi- 
ficance than  simple  geometric  solidity. 

The  relation  between  values  of  CT^and  geometric  solidity 
is  directly  influenced  by  blade  thickness.  Very  thin  blades 
have  geometric  solidities  almost  exactly  equal  to  G" .  With 
thicker  blades,  the  geometric  solidity  is  always  less.   The 
reason  for  this  may  be  illustrated  by  examples  of  two  extremes. 
If  a  contour  in  the  near  circle  plane  is  transformed  confor- 
mally into  a  flat  plate,  the  chord  length  of  the  flat  plate 
will  be  nearly  twice  the  diameter  of  the  original  contour. 
However,  if  no  transformation  is  made  at  all,  an  "airfoil"  of 
100$  thickness  results,  whose  chord  length  is  equal  to  the 
diameter.  Both  cascades  would  have  the  same  value  of  (T* ,  but 
the  geometric  solidity  of  the  latter  case  would  be  half  that 
of  the  former.  Between  these  two  extremes,  finite  thickness 
airfoils  in  cascade  have  solidities  less  than  <r"  . 

2.4  The  Turning  Parameter 

The  turning  parameter,  B0,  is  a  measure  of  the  change  in 


direction  of  the  flow  passing  through  the  blade  row.  It  is 
actually  a  nondimensional  form  of  the  downwash  velocity  induced 
by  the  blades.  There  is  a  very  simple  relation  between  B0  and 
the  circulation  around  each  blade.  In  fact,  for  a  constant 
mean  flow  angle,  this  relation  is  a  simple  direct  proportion. 
This  convenient  state  of  affairs  is  the  reason  why  the  two 
flow  parameters  were  expressed  in  their  present  form. 

In  isolated  airfoil  theory,  the  method  of  introducing  lift 
or  circulation  on  an  airfoil,  is  to  add  to  the  complex  poten- 
tial, a  terra  representing  a  vortex  of  suitable  strength.  How- 
ever, mere  addition  of  a  vortex  will  produce  a  distortion  in 
the  contours  in  the  near  circle  plane,  since  these  contours 
are  not  exactly  circles.  In  addition,  an  infinite  stack  of 
vortices  is  necessary  to  produce  circulation  around  each  of  the 
blades  in  the  cascade.  To  avoid  these  resulting  undesirable 
perturbations  in  what  will  ultimately  become  the  airfoil  shape, 
and  indirect  method  of  adding  circulation  is  used. 

In  this  method,  a  transformation  is  made  from  the  near 
circle  plane  to  the  circle  plane.  An  equation  is  developed 
which  transforms  the  flow  in  the  near  circle  plane  into  the 
flow  around  a  perfect  circle,  centered  on  the  origin.  The  cir- 
culation term  may  then  be  introduced  without  fear  of  distor- 
tion. The  transformation  is  then  made  back  to  the  near  circle 
plane.   Since  the  same  equation  is  used  to  transform  both  ways, 
the  original  flow  contour  remains  unchanged;  however,  the 
velocities  at  each  point  are  now  different,  due  to  the  circu- 

10 


lation. 

2.5  The  Blade  Shape  Parameters 

The  final  transformation  is  from  the  near  circle  plane 
to  the  cascade  plane.  This  is  done  by  the  Joukowski  method 
of  shifting  the  axes  slightly  in  such  a  way  that  the  flow  con- 
tour transforms  into  an  airfoil  shape .  Two  arbitrary  para- 
meters were  chosen  to  control  the  blade  shape.  One  of  these 
is  called  the  thickness  parameter,  and  is  given  the  symbol  £,  . 
This  quantity  is  generally  related  to  the  thickness,  and  also 
to  nose  curvature,  but  not  in  any  fixed  relationship.  Large 
values  of  £,    produce  thick  airfoils,  and  small  values  produce 
thin  airfoils.  Figure  2  is  an  illustration  of  several  airfoils, 
showing  the  effect  of  C*   .  The  limiting  value  of  0  for  £  does 
not  however  produce  an  airfoil  of  zero  thickness.  Due  to  the 
slight  distortion  of  the  contours  in  the  near  circle  plane, 

6*  -  0   may  produce  the  mathematically  significant,  but  physi- 
cally useless  case  where  the  "lower"  surface  of  the  airfoil 
crosses  through  the  "upper"  surface.  £  then,  can  be  small  but 
not  zero,  and  the  airfoils  produced  can  be  thin  but  not  of  zero 
thickness. 

The  other  parameter  is  the  shape  or  symmetry  parameter  K. 
The  effect  of  variations  in  K  may  be  seen  in  Figure  3.  Gener- 
ally, K^O  produces  a  blunt  leading  edge,  and  cusped  trailing 
edge  similar  to  a  Joukowski  airfoil.  K=l  will  produce  a  sym- 
metrical shape  which  is  generally  a  cambered  ellipse.  As  K 
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increases  from  zero  to  one,  the  point  of  maximum  thickness 
moves  from  the  quarter  chord  to  the  half  chord  point,  and  the 
trailing  edge  rounds  off »     Values  of  K  larger  than  one  will 
produce  airfoils  facing  the  wrong  way  (Trailing  edge  into  the 
flow)  and  hence  must  not  be  used. 

2.6  Stagnation  Points 

During  the  transformations  made  in  this  analysis,  the 
stagnation  points  serve  as  reference  points  to  locate  the  re- 
lation of  the  contours  in  the  various  plane So   For  this  reason 
there  is  always  a  stagnation  point  at  the  leading  and  the 
trailing  edge*     Having  a  stagnation  point  at  the  leading  edge 
is  sometimes  called  the  Theodorsen  condition .  This  means  that 
in  all  the  cascades  developed  by  this  method,  the  airfoils  are 
operating  at  the  ideal  angle  of  attack o  Variations  in  lift 
then,  are  achieved  by  variations  in  the  camber  of  the  blades. 
Figure  1  graphically  illustrates  this  facto  Two  typical 
cascades  are  drawn,  with  the  principal  variation  being  in  turn- 
ing parameter,  BQ.  The  values  of  £  were  adjusted  to  give 
the  same  blade  thickness »   The  difference  in  camber  may  be 
clearly  seen*. 

The  condition  of  a  stagnation  point  at  the  trailing  edge 
is  called  the  Kutta  condition.  This  fact  has  particular  sig- 
nificance when  considering  the  blade  shapes  where  K=lo  As 
pointed  out  previously,  these  shapes  are  generally  cambered 
ellipses.  It  is  a  well  known  fact  that  flow  around  an  ellipse 
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of  a  given  thickness  ratio  produces  a  lower  velocity  peak  than 
the  flow  around  any  other  shape  of  the  same  thickness  ratio. 
This  class  of  blades  then  should  exhibit  the  lowest  loading 
attainable  for  any  possible  fairing  shape.  Of  course,  any 
attempt  to  duplicate  this  fairing  shape  in  an  actual  cascade 
is  doomed  to  failure  because  of  the  rounded  trailing  edge. 
However,  this  case  still  has  significance  in  two  respects. 
First,  the  performance  of  a  cambered  ellipse,  with  the  Kutta 
condition  mathematically  imposed,  represents  a  useful  limiting 
case  of  performance  due  to  fairing  shape.  Secondly,  a  blade 
similar  to  this,  with  modifications  to  produce  a  sharp  trail- 
ing edge,  would  have  very  nearly  the  same  performance.   The 
critical  area  of  flow  near  the  leading  edge  would  not  be  af- 
fected much  by  modifications  at  the  rear  of  the  airfoil. 

Figure  4  is  a  plot  of  the  pressure  distribution  around 
the  airfoils  of  Figure  3»   The  above  premise  is  verified,  as 
the  negative  pressure  peak  can  be  seen  to  decrease  as  the 
value  of  K  is  increased.  Note  that  the  pressure  distribution 
aft  of  the  peak  is  essentially  linear  to  the  trailing  edge. 
Reference  3,  p»  354  states  that  this  type  of  pressure  distri- 
bution is  typical  of  decelerating  cascades.  Note  also  that 
the  slope  of  this  line  becomes  less  as  K  is  increased,  indi- 
cating a  less  adverse  pressure  gradient  in  the  boundary  layer. 
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3.0  Prediction  of  Separation 

3.1  Pressure  Gradient 

The  phenomenon  which  ultimately  limits  the  attainable 
pressure  rise  of  a  cascade  is  separation.  Flow  separation 
occurs  when  flow  in  a  boundary  layer  experiences  an  adverse 
(positive)  pressure  gradient.  This  adverse  gradient  tends  to 
decrease  the  velocities  in  the  boundary  layer.  Ultimately, 
the  point  is  reached  where  the  velocities  close  to  the  surface 
have  been  reduced  to  zero.  At  this  point,  flow  separation 
occurs.  Beyond  this  point,  some  velocities  become  negative, 
and  a  region  of  reverse  flow  exists.  The  resultant  effects 
cause  high  drag  and  loss  of  circulation  around  each  blade 
element.  Many  experiments  have  verified  that  any  adverse 
pressure  gradient,  however  small,  will  cause  separation  if  it 
acts  long  enough.   The  concept  of  a  gradient  and  the  distance 
over  which  it  acts,  leads  one  naturally  to  suspect  pressure 
difference  as  a  significant  parameter  in  predicting  separation, 
The  gradient,  integrated  over  the  distance  is  mathematically 
equal  to  the  difference  in  pressure .  It  was  previously  shown, 
furthermore,  that  the  pressure  distribution  over  typical  air- 
foils in  this  system  is  approximately  linear  aft  of  the  nega- 
tive pressure  peak. 

3.2  Velocity  Considerations 

From  another  point  of  view,  one  might  examine  the  velo- 
cities of  flow  through  a  cascade .   From  the  inlet  velocity 
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described  previously,  velocity  is  Drought  to  zero  at  the 
front  stagnation  point,  accelerated  to  a  local  maximum  some- 
where on  the  surface  of  the  blades,  reduced  toward  zero  at 
the  trailing  edge,  and  finally  settles  down  to  the  outlet 
velocity  far  behind  the  cascade o  Besides  the  inlet  and  outlet 
velocities,  only  one  other  velocity  is  really  significant,  and 
that  is  the  local  maximum,  Vmax  .   From  these  three  quantities, 
only  two  independent  diraensionless  parameters  may  be  construc- 
ted. If  the  inlet  velocity  and  Vmax  were  used,  this  parameter 
would  be  related  to  the  flow  around  the  forward  part  of  the 
airfoil.  Since  this  flow  involves  a  negative  pressure  grad- 
ient, it  is  probably  not  very  significant  with  respect  to 
separation.  Similarly,  the  quantities  V    and  the  outlet 
velocity,  called  V2»  are  related  to  the  flow  around  the  rear 
part  of  the  airfoil,  which  involves  a  positive  pressure  grad- 
ient. These  two  parameters  then,  combined  in  some  diraension- 
less manner,  should  act  as  an  indicator  for  the  tendency 
toward  separation. 

3-3  Diffusion  Factors 

From  its  many  experimental  studies  of  flow  through  cas- 
cades, the  NASA  has  arrived  at  essentially  the  same  conclusion. 
They  have  used  several  diraensionless  combinations  of  Vmax  and 
V2  to  define  blade  loading  limits.  One  of  these  is  the 
"Local  Diffusion  Factor"  which  is: 
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This  factor  is  described  in  Reference  1  and  supporting  data 
indicates  that  separation  is  likely  to  occur  for  values  of  D^ 
greater  than  0,5 « 

Since  the  velocity  distribution  around  the  blades  must 
be  known  to  find  values  of  D^  ,  it  is  necessary  to  use  experi- 
mental data,  or  an  analytical  treatment  such  as  the  one  des- 
cribed in  this  investigation.  To  get  around  this  difficulty, 
the  NASA  also  defines  an  approximation,  D,  which  may  be  ob- 
tained from  more  easily  determined  quantities .  D,  is  called 
the  "Diffusion  Factor"  and  is  given  by: 

D  -  1-  Cos£l  +  Cos/Ql  (Tan/21  -  Tan/02) 
Cos  /%     2  ys 

It  is  seen  that  this  quantity  can  be  determined  from  the  inlet 
and  outlet  velocities,  and  the  solidity  of  the  blade  row,  The 
supporting  data  indicates  that  separation  is  likely  to  occur 
for  values  of  D  greater  than  0o6„ 

Since  the  methods  used  in  this  investigation  furnish  all 
of  the  flow  quantities  described  above,  as  well  as  the  com- 
plete velocity  distribution  around  each  blade,  the  opportunity 
to  compare  the  two  NASA  parameters  is  presented.   As  will  be 
described  later,  the  correlation  between  these  two  parameters 
is  not  too  good  for  the  solidities  considered  hereo 
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3«4  The  Boundary  Layer  Loading  Parameter 

From  the  above  development,  it  may  be  seen  that  the  para- 
meter to  indicate  separation  should  involve  the  pressure  ratio 
across  the  rear  part  of  the  airfoil,  and  be  some  dimensionless 
combination  of  the  quantities  V2  and  Vfflaxo   The  NASA  has  in- 
dicated that  most  compressors  operate  in  the  range  where 
Reynolds  number  effects  are  insignificant „  Therefore,  it  is 
felt  that  no  particular  inaccuracies  result  from  omitting 
consideration  of  Reynolds  number . 

In  order  to  avoid  using  one  parameter  for  pressure  dis- 
tribution, and  another  for  the  tendency  toward  separation, 
the  quantity  S  was  selected  where; 

s  -  (  W  v2>2 

This  quantity  is  referred  to  as  the  boundary  layer  loading 
parameter,  and  is  quite  versatile  in  application .   It  is 
directly  related  to  the  pressure  coefficient,  and  in  fact  only 
differs  from  that  quantity  by  unity »  S  also  has  a  constant 
functional  relationship  with  the  NASA  parameter  D^o   The  value 
S  =  4*25  corresponds  to  the  value  of  D^  of  0o5,  and  has  been 
selected  as  the  limit  loading  in  the  boundary  layer  for  blades 
in  this  investigation o 
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4.0  Results  and  Discussion 

Theory  and  experiment  both  point  toward  increasing 
solidity  as  the  method  of  obtaining  increased  pressure  rise 
per  stage •  As  solidity  is  increased,  the  lift  coefficient  on 
each  blade  is  decreased,  and  hence  the  loading  on  the  boundary 
layer  is  decreased  alsoo  However,  as  the  solidity  is  increased, 
there  is  an  adverse  effect  due  to  the  blockage  caused  by  the 
finite  thickness  of  the  blades  plus  the  boundary  layer » 
Finally,  the  point  is  reached  where  the  beneficial  effects  of 
decreased  blade  loading  are  just  counterbalanced  by  the  ad- 
verse effects  of  blockage o  The  quantitive  determination  of 
this  point  was  one  of  the  major  objectives  of  this  investiga- 
tion. 

4.1  Organization  and  Scope 

Since  actual  cascades  are  designed  on  the  basis  of  blade 
thickness  and  geometric  solidity,  it  was  felt  that  the  results 
of  this  investigation  ought  to  be  organized  on  this  basis, 
rather  than  on  the  basis  of  the  mathematical  parameters  £•  and 
O^   .  The  value  of  blade  thickness  is  not  directly  predictable 
from  the  input  parameter  £  •  Therefore,  the  basic  computer 
program  was  modified  to  use  an  iterative  scheme  to  achieve  the 
desired  values  of  blade  thickness •  Details  are  given  in 
Appendix  B,  which  describes  the  basic  computer  program  and  the 
modifications .  Seven  blade  thicknesses  were  selected  to  cover 
the  field  of  possible  blade s0  These  ranged  from  the  limiting 
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case  of  extremely  thin  blades  to  2fy$   thickness  blades  in  four 
percent  increments . 

It  was  decided  to  conduct  this  investigation  using  only- 
one  value  of  average  flow  angle  /So  This  value  was  chosen  to 
be  45°,  which  was  a  somewhat  arbitrary  decision,  but  not  with- 
out some  sound  theoretical  basis.  It  is  stated  in  Reference 
4  that  maximum  efficiency  may  be  obtained  when  the  average 
flow  angle  is  near  45  °»  and  many  compressors  are  built  accord- 
ingly. 

In  spite  of  reducing  the  parameters  to  four  by  the  above 
decision,  over  1470  separate  cascades  were  generated  and  ana- 
lyzed for  this  investigation.  Additional  parameters  could 
have  been  introduced  to  give  more  control  over  the  blade  shape, 
but  they  would  have  greatly  magnified  the  problem.  As  the 
data  will  show,  additional  blade  shape  parameters  would  not 
have  affected  the  significance  of  the  results  very  much. 

4.2  Series  1  Plots 

The  first  plots  of  data  from  the  computer  run  were  made 
of  boundary  layer  loading  vs.  solidity,  for  airfoils  of  con- 
stant blade  thickness.   Five  values  of  the  shape  parameter  K 
were  used,  and  data  was  obtained  at  six  values  of  turning 
parameter,  BQ0   This  meant  that  there  were  30  graphs  in  series 
one.   Figures  6  through  9  are  representative  graphs  from  this 
series. 

It  was  noted  with  satisfaction  that  the  curves  showed  a 
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pronounced  minimum  loading  as  solidity  was  increased.  The  one 
exception  was  the  case  of  the  extremely  thin  airfoils,  for 
which  the  curve  appeared  to  be  asymptotic  to  some  low  value  of 
loading.  This  is  not  unreasonable,  since  the  blockage  effect 
would  be  very  small  for  this  case.  The  case  of  thin  airfoils 
then,  represent? a  sort  of  limiting  value  for  cascades.  The 
point  at  which  minimum  loading  occurred  for  a  given  thickness 
airfoil  was  called  the  optimum  solidity.. 

It  may  be  noted  that  while  the  value  for  loading  at  an 
optimum  point  could  be  determined  with  accuracy,  the  corres- 
ponding solidity  for  that  point  was  subject  to  some  uncertainty. 
Therefore,  another  point  was  introduced  where  the  solidity 
could  be  determined  with  greater  certainty.  This  was  called 
the  design  point,  and  is  defined  as  the  point  toward  lower 
solidity  where  the  loading  equals  1.10  times  the  minimum  value. 
There  is  further  justification  in  this  step  from  the  fact  that 
the  basic  analysis  neglects  boundary  layer.  The  blockage 
effect  of  a  blade  plus  boundary  layer  must  be  greater  than  the 
effect  due  to  the  blade  alone.  Therefore,  the  true  optimum 
solidity  must  be  somewhat  less  than  the  apparent  optimum  based 
on  non-viscous  flow.  The  exact  relation  cannot  be  found  with- 
out a  detailed  analysis  of  the  boundary  layer.  This  simple 
method  sidesteps  this  formidable  problem,  and  is  at  least 
qualitatively  in  the  right  direction. 

For  the  case  of  higher  values  of  turning  parameter  and 
thick  blades,  there  was  some  scatter  in  the  data.   Since  the 
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data  is  all  analytical,  and  should  therefore  be  reasonably- 
regular,  some  effort  was  made  to  find  the  source  of  the  scat- 
ter. The  relationship  between  the  parameters  C  and  C/S  was 
suspected  first.  Comparative  plots  of  these  two  parameters 
however,  showed  that  they  varied  in  a  regular  fashion  with 
blade  thickness.  It  was  then  decided  that  the  scatter  was  due 
to  the  fact  that  velocities  were  measured  only  at  40  points 
along  the  airfoil.  If  the  true  maximum  velocity  occurred  be- 
tween these  finite  points,  then  the  value  for  boundary  layer 
loading  parameter  would  be  in  error.  Remembering  the  fact 
that  this  error  could  be  only  on  the  low  side  was  an  aid  in 
fairing  the  curves  through  the  scattered  points. 

It  was  noted  that  the  data  became  more  regular  as  the 
value  of  K  approached  one.  For  these  airfoils,  either  the 
flow  was  more  regular  about  the  nose,  with  no  sharp  peaks  in 
velocity,  or  the  point  of  maximum  velocity  happened  to  coin- 
cide more  often  with  one  of  the  finite  points.  Examination 
of  the  pressure  plots  of  Figure  4  points  toward  the  former 
as  the  cause.  With  high  values  of  K,  the  pressure  changes  are 
seen  to  be  more  gradual,  and  of  lower  magnitude.  Because  of 
this,  the  actual  maximum  velocity  was  always  close  to  the 
velocity  recorded  at  one  of  the  40  points • 

Figure  10  is  an  illustration  of  the  effect  of  K  in  reduc- 
ing the  boundary  layer  loading  for  a  given  value  of  turning 
parameter.  With  all  other  factors  equal,  the  loading  steadily 
decreased  as  the  value  of  K  increased  toward  unity=   Note  that 
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the  effect  of  K  was  much  greater  lor  the  thick  blades.   For 
the  case  of  extremely  thin  blades,  K  had  no  effect  at  all, 
which  was  not  an  unexpected  result . 

Jm3  Series  2  Plots 

Since  the  original  series  involved  30  graphs,  some  effort 
was  made  to  summarize  this  data.   The  result  of  this  effort 
was  Figures  11  through  15,  which  collect  all  the  data  for  a 
particular  value  of  K  on  a  single  page.  Boundary  layer  load- 
ing parameter  and  solidity  are  plotted  along  the  coordinates 
as  before,  but  now  the  solidity  is  design  solidity.  In  effect, 
all  the  design  points  for  a  particular  value  of  K  have  been 
plotted.   Points  of  constant  values  of  turning  parameter  have 
been  connected,  as  have  points  of  constant  airfoil  thickness. 
The  result  is  that  four  interdependent  parameters  are  shown 
for  each  point  on  the  chart.   Any  chosen  point  then,  repre- 
sents a  particular  value  of  blade  thickness  and  turning  para- 
meter, with  the  corresponding  boundary  layer  loading  and 
design  solidity  for  the  cascade. 

Figures  16  through  20  are  plots  made  in  exactly  the  same 
manner,  except  that  optimum  points  were  used  instead  of  design 
points.   It  was  necessary  to  change  the  scale  for  solidity 
along  the  horizontal  axis,  since  optimum  points  represent 
higher  values  of  solidity  in  every  case.   In  addition  some 
difficulty  was  experienced  with  irregularity  in  the  data. 
This  was  due  at  least  in  part  to  the  uncertainty  in  determin- 
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ing  the  value  of  solidity  at  the  optimum  point*  This  optimum 
point  data  was  included  as  a  limiting  case,  to  illustrate  the 
maximum  solidities  to  be  used  when  the  boundary  layer  is  neg- 
ligibly thin. 

The  choice  of  solidity  as  a  coordinate  is  slightly  mis- 
leading.  It  might  be  inferred  from  this  choice  that  solidity 
should  be  used  as  a  design  input •   However,  solidity  was 
plotted  in  this  manner  purely  as  a  convenience  in  transferring 
data  from  the  series  one  graphs.  In  design  work,  the  two 
"inputs"  could  be  boundary  layer  loading,  and  turning  para- 
meter.  The  "output"  information  would  then  be  blade  thickness 
and  solidity.  An  attempt  was  made  to  plot  boundary  layer 
loading  parameter  vs.  turning  parameter,  but  the  resulting 
blade  thickness  and  solidity  lines  were  so  nearly  parallel 
that  the  information  was  difficult  to  read. 

4.4  Comparison  of  NASA  Separation  Parameters 

The  data  on  the  series  2  plots  contains  all  of  the  infor- 
mation necessary  to  compute  both  NASA  separation  parameters. 
This  affords  an  opportunity  to  compare  these  two  parameters 
directly.   The  boundary  layer  loading  parameter  S,  and  the 
NASA  parameter  D-^  have  a  constant  relationship  such  that: 

B1   =  1  -  \_ 

Thus  a  D-i  equal  to  0.5  corresponds  to  a  value  of  S  of  4»25. 
This  value  is  entered  as  a  heavy  line  on  all  series  2  plots. 
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The  approximation  b,  is  related  to  the  parameters  B,  BQ, 
and  c/S  by  the  following  relationship; 


D.  ,-lh  +  faft-  if  * 

f  I  +  (tan/3  r  %-f  *  SKf^/S^ff 

Thus,  since  (3   has  been  fixed  at  one  value,  each  value  of  B0 
and  c/S  will  result  in  a  particular  value  for  Do   This  expres- 
sion for  D  shows  that  it  cannot  be  a  good  approximation  over 
a  wide  range  of  solidities,  since  it  indicates  that  loading 
will  always  decrease  whenever  solidity  is  increased .  It  has 
already  been  demonstrated  that  as  solidity  increases  beyond 
the  optimum  point,  loading  increases  also  due  to  the  blockage 
effect. 

Lines  corresponding  to  a  value  of  D=.5  and  D*=.6  have  been 
plotted  on  Figures  11  through  20  as  dashed  lines<>   A  note  of 
caution  is  necessary  in  using  the  data  on  these  figures.  It 
must  be  remembered  that  there  is  really  only  one  solidity 
shown  on  each  chart,  and  that  is  the  design  (or  optimum)  soli- 
dity.  The  fact  that  solidity  is  shown  as  a  variable  is  due 
solely  to  the  fact  that  each  combination  of  blade  thickness 
and  turning  parameter  demands  a  particular  value  of  solidity 
at  the  design  (or  optimum)  point. 

If  the  parameter  D  is  a  true  approximation  for  D^.  then 
the  lines  corresponding  to  D=.5  and  D^-c5  should  be  nearly 
the  same.   Examination  of  the  figures  shows  that  there  is  a 

24 


fair  match  for  the  case  of  optimum  solidity.   Generally,  D 
gives  a  value  of  loading  which  is  nigher  than  D-p  but  the 
error  is  only  around  10$,  which  is  not  unreasonable.  For  the 
case  of  design  solidity,  there  is  a  greater  discrepancy.  For 
very  thick  blades,  D  still  indicates  a  value  of  loading  higher 
than  D^>  but  this  relation  reverses  so  that  for  thin  blades 
D^  indicates  the  higher  loading , 

The  data  in  Reference  1  states  that  separation  is  likely 
to  occur  for  values  of  D^  greater  than  ,5  or  for  values  of  D 
greater  than  .6,  This  in  itself  is  an  indication  of  a  mis- 
match between  D  and  D^.  If  the  occurrence  of  separation  is 
taken  as  the  matching  criterion,  then  the  line  corresponding 
to  D=.6  should  match  the  Dt  =  ,5  line.   For  the  case  of  optimum 
solidity  it  may  be  seen  that  the  matching  is  very  poor,  where 
D  indicates  a  much  higher  loading  than  D-,,   At  design  solidity 
the  matching  is  better,  but  D  still  indicates  values  of  load- 
ing which  are  around  18%   higher  than  D^,  Extrapolating  this 
trend  leads  one  to  suspect  that,  for  some  value  of  solidity 
less  than  design  solidity,  the  two  parameters  D=,6  and  D^=,5 
would  probably  match  quite  closely.   From  this,  it  is  not  un- 
reasonable to  expect  that  the  NASA  obtained  its  experimental 
data  from  cascades  whose  solidity  was  less  than  the  design  or 
optimum  as  defined  in  this  report,   A  few  checks  of  source 
data  for  Reference  1  have  verified  this  premise.   Solidities 
were  generally  1.0  for  cascades  with  blade  thicknesses  of 
eight  and  ten  percent,   A  check  of  the  series  2  plots  will 
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reveal  that  this  is  less  than  the  aesign  or  optimum  point  for 
any  fairing  shape . 

4»5  Comparison  of  Blade  Thickness  and  Fairing  Shape  Effects 
Figures  11  through  20  clearly  show  the  gains  that  can  be 
made  in  turning  parameter,  at  no  increase  in  loading,  by 
using  thinner  blades  at  the  corresponding  higher  solidities. 
However,  it  was  desired  to  compare  the  effects  of  changes  in 
fairing  shape  and  blade  thickness  on  the  same  chart.   Figure 
21  was  made  using  the  design  point  data,  and  Figure  22  was 
made  with  data  from  the  optimum  points •   The  value  of  S  ■  4*25 
was  selected  as  the  best  estimate  of  maximum  permissible  blade 
loading.   The  values  plotted  were  turning  parameter  vs.  solid- 
ity at  this  constant  value  of  blade  loading.   This  resulted  in 
lines  which  represented  varying  blade  thickness  at  a  constant 
value  of  K.  Lines  of  constant  blade  thickness  were  then  drawn 
in.   It  is  interesting  to  note  that  these  lines  are  nearly 
vertical.  This  means  that  the  value  of  design  or  optimum 
solidity  for  a  particular  blade  thickness  is  relatively  inde- 
pendent of  blade  shape.   This  result  was  not  directly  predic- 
table from  theory,  and  may  prove  useful  in  planning  experi- 
mental investigations. 

As  blade  thickness  is  decreased,  it  is  clearly  seen  that 
the  importance  of  fairing  shape  decreases  correspondingly. 
When  blade  thickness  reaches  about  four  percent,  there  is 
almost  no  discernible  difference  due  to  fairing  shape.   To 
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compare  the  relative  effects  of  fairing  shape  and  thickness 
changes,  consider  the  following  example*  Select  a  blade  thick- 
ness of  12%  with  K  =  .010,  which  is  a  relatively  poor  fairing 
shape  with  a  blunt  nose.  For  this  shape  at  the  proper  solid- 
ity, a  turning  parameter  of  .74  is  possible  without  separation. 
If  K  ■  .990  the  shape  is  nearly  the  best  possible. 

Making  this  change  increases  the  permissible  turning 
parameter  to  .83,  or  a  12%  improvement .  This  same  improvement 
could  be  obtained  with  the  original  shape  by  decreasing  the 
thickness  to  a  little  less  than  3%  and  increasing  the  solidity 
slightly.  If  the  blades  could  be  decreased  to  4%,  with  the 
corresponding  increase  in  solidity,  the  turning  parameter 
could  be  increased  to  .92  without  changing  the  loading.  This 
is  a  24%  improvement,  or  double  the  best  improvement  that 
could  be  obtained  by  fairing  shape  changes  alone.  As  a  limit- 
ing value,  if  the  blades  could  be  made  infinite smally  thin, 
the  increase  in  turning  parameter  could  be  as  high  as  50%. 

It  must  be  remembered  that  these  performance  figures  are 
for  cascades  operating  only  at  the  ideal  angle  of  attack.  It 
may  well  be  that  these  increases  in  peak  performance  are  at- 
tained with  some  sacrifices  in  flexibility „   This  would  be  of 
no  consequence  in  a  constant  output  device,  but  for  general 
application,  the  problem  of  operating  off -de sign  must  be  con- 
sidered.  The  analytical  method  described  in  this  study  can 
be  modified  to  consider  the  case  of  operating  at  angles  of  at- 
tack other  than  the  ideal •   This  matter  is  discussed  more 
fully  in  the  section  on  suggestions  for  further  work. 
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5.0  Conclusions 

It  is  concluded  that  the  method  of  generating  cascades 
as  described  in  this  report  is  a  workable  one,  and  provides 
useful  results.  Cascades  are  generated  as  complete  entities, 
which  eliminates  the  need  for  separating  the  effects  of  the 
various  design  parameters.   Use  of  a  high  speed  digital  com- 
puter means  that  large  numbers  of  cascades  can  be  generated 
and  analyzed  quickly.  This  allows  exploration  over  a  wide 
range  of  design  parameters. 

Solidity  in  experimental  cascades  has  often  been  arbi- 
trary.  This  investigation  has  revealed  that  the  matter  of 
solidity  is  of  paramount  importance.   For  any  combination  of 
fairing  shape,  thickness  and  camber,  there  is  a  particular 
value  of  solidity  which  will  result  in  the  lowest  loading  of 
the  boundary  layer.   This  loading  can  increase  quite  rapidly 
as  the  solidity  is  varied  from  the  proper  value. 

Data  obtained  in  this  investigation  has  yielded  the  un- 
expected result  that  the  optimum  value  of  solidity  referred 
to  above  is  relatively  independent  of  the  fairing  shape  of  the 
blades.   An  important  application  of  this  result  will  be  in 
planning  experimental  investigations  of  cascade  performance. 
The  actual  values  of  solidity  given  in  this  report  may  be  used, 
with  the  probability  that  they  will  be  close  to  the  optimum 
values  for  any  sort  of  physical  blade  shape. 

Three  parameters  intended  as  indicators  of  the  tendency 
toward  flow  separation  were  compared  in  this  report.   As  a 
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result  of  this  comparison,  it  is  concluded  that  the  parameter 
S,  as  defined  in  this  report,  is  a  useful  indicator  in  this 
respect.   It  is  further  concluded  that  the  value  S  =  4.25  is 
an  acceptable  design  limit,  which  is  as  consistent  with  NASA 
data  as  this  data  is  consistent  with  itself. 

With  respect  to  the  problem  of  increasing  the  performance 
of  airfoils  in  cascades,  it  is  concluded  that  the  matter  of 
blade  thickness  is  of  much  greater  importance  than  the  fairing 
shape.   The  example  was  given  of  a  cascade  with  blades  of  12% 
thickness.   Possible  increases  in  turning  parameter  for  this 
cascade  due  to  improvements  in  fairing  shape  were  shown  to  be 
on  the  order  of  12%.   Increases  in  turning  parameter  produced 
by  decreasing  the  blade  thickness  to  4%  with  a  corresponding 
increase  in  solidity,  were  shown  to  be  24%,  or  twice  the  pre- 
vious amount.   The  absolute  limit  of  performance  increase  by 
this  method  would  be  given  by  the  case  of  negligibly  thin 
blaaes.   For  this  case,  the  total  increase  in  turning  parameter 
would  be  50%. 

It  was  pointed  out  in  this  study  that  increasing  the 
performance  of  cascades  by  using  thinner  blades  and  higher 
solidities  would  probably  reduce  the  flexibility  of  performance, 
The  effects  of  operating  off -design  were  not  a  part  of  this 
study  since  the  method  of  generating  the  cascades  produced 
airfoils  operating  only  at  the  ideal  angle  of  attack.   It  was 
stated  that  modifications  could  be  made  in  the  method  to  con- 
sider the  off -design  case,  but  this  was  not  done  due  to  time 
limitations. 
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6.0  Suggestions  for  Future  Work 

This  investigation  represents  only  a  beginning  in  the 
field  of  theoretical  design.  With  the  basic  method  and  com- 
puter programs  available,  it  should  be  a  simple  matter  to  ex- 
tend the  scope  of  investigations  similar  to  this  one.   Inves- 
tigations ought  to  be  conducted  for  a  number  of  average  flow 
angles,  so  as  to  cover  the  whole  range  of  possibilities  in 
cascade  design.  Reference  2  for  example,  has  some  work  on  an 
average  flow  angle  of  zero,  which  represents  the  special  case 
of  an  impulse  turbine.   There  is  no  reason  to  prevent  con- 
sidering negative  average  flow  angles  either.   This  would  be 
the  case  of  turbines  rather  than  compressors,  but  the  same 
method  and  computer  program  would  apply. 

One  might  suppose  that  since  the  flow  is  accelerated 
through  a  turbine,  the  boundary  layer  experiences  a  negative 
pressure  gradient  and  separation  is  not  a  problem.  This  is 
not  exactly  true,  however.  There  are  actually  two  pressure 
differences  to  consider  in  flow  through  a  cascade.  One  is  the 
effect  perpendicular  to  the  stream,  which  is  a  result  of  the 
lift  being  developed  by  each  blade.  The  other  is  the  longi- 
tudinal effect  parallel  to  the  stream  which  is  the  result  of 
differences  in  pressure  from  the  front  to  the  rear  of  the 
blade  row.   Flow  separation  depends  on  the  combination  of  these 
two  effects.  Thus,  in  decelerating  cascades,  the  longitudinal 
effect  reduces  the  allowable  lift  from  each  blade.   In  tur- 
bines, the  longitudinal  effect  increases  the  allowable  lift. 
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If  the  lift  is  increased  enough  to  overcome  the  aid  of  the 
longitudinal  effect,  flow  separation  will  still  occur.  A 
series  of  studies  of  this  matter  could  establish  the  exact 
relation  between  these  two  gradients  and  the  occurrence  of 
separation. 

The  analytical  methods  of  this  report  would  be  a  help  in 
studying  the  occurrence  of  separation  through  solution  of  the 
boundary  layer  equations.  The  detailed  pressure  distributions 
provided  by  the  computer  at  least  open  an  avenue  toward  solu- 
tion of  the  boundary  layer  equations  by  direct  integration. 

In  any  future  work,  it  would  be  convenient  to  consider 
the  velocity  at  more  points  near  the  leading  edge  of  the  blades. 
This  would  allow  more  accurate  determination  of  the  actual 
maximum  velocity,  and  should  eliminate  some  of  the  scatter 
found  in  the  data  for  this  investigation.  The  necessary  modi- 
fications to  the  computer  program  are  not  difficult,  and  only 
lack  of  time  prevented  their  incorporation  in  this  investiga- 
tion. 

Another  area  of  investigation  would  be  a  change  in  the 
basic  transformation  which  would  allow  a  finite  angle  at  the 
trailing  edge  of  the  blades.   This  would  make  the  transforma- 
tions similar  to  the  Karmann-Trefftz  transformations  then, 
rather  than  the  Joukowski.  Eliminating  the  cusp  or  radius  at 
the  trailing  edge  then,  would  perhaps  yield  blade  shapes  which 
would  be  more  physically  useful.   As  pointed  out  earlier,  some 
of  the  blade  shapes  developed  in  this  investigation,  while 
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mathematically  useful,  were  physically  poor  because  of  the 
rounded  trailing  edge.   These  improved  transformations  would 
be  particularly  useful  in  the  case  of  symmetrical  impulse 
blading,  where  a  finite  angle  is  desired  at  both  the  leading 
and  the  trailing  edge. 

As  mentioned  previously,  it  would  be  interesting  and  per- 
haps useful  to  consider  the  off -design  case,  where  the  blades 
are  operating  at  some  angle  of  attack  other  than  the  ideal. 
In  principle,  the  method  would  be  to  generate  the  airfoils, 
then  go  back  to  the  original  complex  potential  and  change  the 
average  flow  angle.   Keeping  the  same  transformations  would 
preserve  the  blade  shape,  but  the  circulation  would  have  to 
be  adjusted  to  keep  the  Kutta  condition.   The  idea  seems  simple 
enough,  but  no  doubt  there  will  be  many  difficulties  with  the 
details. 

Still  further  work  could  be  done  in  investigating  other 
shapes  for  the  camber  line  of  the  biadeso   Generally,  the 
camber  line  for  the  blades  developed  in  this  investigation  was 
the  slightly  flattened  arc  of  a  circle .   The  generation  of 
other  camber  lines  involves  the  expression  for  the  flow  poten- 
tial in  the  near  circle  plane .   In  this  plane,  a  single  doub- 
let generates  the  flow  contour  which  transforms  into  one  of 
the  blades  of  the  cascade.   If  instead,  a  series  of  doublets 
were  used;  by  varying  the  strength  and  distribution  of  these 
doublets,  great  control  could  be  exercised  over  the  camber 
line  and  blade  shape  of  the  airfoils .   The  ultimate  in  this 
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direction  would  be  where  the  input  is  a  desired  pressure 
distribution,  with  the  proper  camber  and  blade  shape  being 
automatically  produced.  The  resultant  family  of  cascades 
would  be  as  important  a  milestone  to  cascade  design  as  the 
NASA  6-series  airfoils  were  to  isolated  airfoil  design. 
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Boundary  Layer  Loading  Parameter  vs.  Solidity 
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Boundary  Layer  Loading  Parameter  vs.  Lesign  Solidity 
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Boundary  Layer  Loading  Parameter  vs.  optimum  Solidity 
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Boundary  Layer  Loading  Parameter  vs.  Optimum  Solidity 
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Boundary  Layer  uoading  parameter  vs.  optimum  Solidity 
jj'or  Various  Cascades     Beta  =  45     K  =  .745 


53 


Series  Two 


9.0 


8.0- 


>_V 


*  7.0- 


► 
a 


<£  6.0 


^ 


4  5.o4 


-i 


4.0-:, 


0  "5  0 


2.0 


1.0 


..  . 

B0-1.25N 

24*° 
X      20°/o 

Bo=lW 

/x  16 

8% 

Jj  =  _.6 

V--I5 

47° 

B  =.30 
o     \ 

_      /     ___  ^^ 

jfc.   —     —       / 

D.=  .-"    " 

__5__ 

l 

b°n// 

e0=.5o/^ 

______________ 



.-_-________. 

1.0 


1.5 


2.0      2.5       3.0 


3.5 


4.0 


Optimum  Solidity 

figure  20. 
Boundary  Layer  Loading  Parameter  vs.  uptimum  Solidity 
.'or  Various  Cascades    Beta  =  45     K  =  .990 


54 


l.lO 


1.05" 


1.2     1.4     1.6 

.Design  Solidity 


figure  21. 

Turning  i-arameter  vs.   Design  solidity 
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(   Corresponds   to  NASA  Local  uif fusion  factor  of   .5    ) 
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Turning  Parameter  vs  optimum  oolidity 

ivor  Boundary  Layer  Loading  parameter  of  4.2b 

(  Corresponds  to  NASA  Local  Lif fusion  factor  of  .5  ) 

56 


MATHEMATICAL  ANALYSIS  Ui?1  Xfici  METHuJJ  ui?  DBSIGJli 

This  appenaix  describes  the  general  mathematical 
development  of  the  method  of  designing  airfoils  in  cascale. 
bome  of  the  fine  details,  particularly  on  iterative  proceed- 
ures,  have  been  omitted,   nowever,  the  entire  calculation 
proceedure  is  given  in  the  Tu.uT.tuu*  program  in  the  next  appen- 
dix.  This,  along  with  the  explanations  given  in  the  table 
of  symbols,  should  allow  analysis  of  the  method  to  whatever 
extent  is   necessary. 

Basic  relations  in  the  Near  Circle  rlane 

The  complex  potential  for  the  flow  around  the  infinite 
row  of  closed  bodies  in  the  near  circle  plane  may  be  written 
in  the  form: 

w  =  <f  +  if  =  e~z   +  (A  +  IB)  2Az)  Ala 


With  the  corresponding  complex  velocity: 

dw  =  v  -  iv  =  e~   +  (A  +  iB)  F„(z) 
dz    x     y 


Alb 


where 

Fx(z)  =  iL^XjV)  +  iJ^Ujy)  =  i   coth  z  *2a 

F2(z)  =  K2(x,y)  +  iJ2(x,y)  =  f^  F±(z) 


=   jL  ^   Coth  f  j     =  -  i   csch2  f 


A2b 


By  expanding  these  functions  and  separating  tnem  into 
their  real  and  imaginary  parts,  it  is  found  that: 
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JI(x^  =  -  *  (cosh  f-  Ls  y  )  A3b 


i  te! 

"    rco: 


JoU.y)  =  +  •   y^g  *  sOa  y  )  A5d 

2  *    (cosh  x  -   cos  yJZ  ^ 


ihe  expressions  for  velocity  potential,  stream  function, 
anu.  velocity  components  now  reduce  to: 


x  cos/2  +  y  slnfi     +  AK^Xjy)  -  BJx(x,y)    A4a 
y  =  -  x  sin/S  +  y  cos(3  +  AJ^x.y)  +  BR-(x,y)   A4b 

cos/3  +  Aiv2(x,y)  -  BJ2(x,y)   A4c 


v  = 
x 


v  =  sin/S  -  AJ2(x,y)  -  BR2(xfy)   A4d 


ihese  relations  represent  the  flow  without  circulation 
about  an  infinite  series  of  nearly  circular  bodies  spaced 
at  equal  intervals  S  =  2  It    along  the  y  axis,   ihe  contour  of 
a  typical  body  is  shown  in  .figure  !?. 

Location  of  Stagnation  Points 

.Due  to  the  interference  caused  by  thu  presence  of  the 
infinite  stack  of  bodies,  the  stagnation  points  on  any  body 
are  displaced  slightly  from  the  position  they  would  have 
occupied  if  the  body  were  isolated  in  the  flow,   r'or  this 
reason,  some  computation  is  necessary  to  locate  these  points. 


designating  the  front  ana  rear  stagnation  points  by  the 

subscripts  A  and  B  respectively,  ana  noting  the  type  of  polar 

symmetry  involved,  the  following  facts  are  observed.   At  the 

two  stagnation  points,  both  velocity  components  vanish  by 

definition,   ^.lso,  by  reason  of  the  odd  symmetry  involved, 

x^  =  -  xB  ,   yA  =  -  y^  ,   and  fA   =-^   .   however,  since 

the  points  A  and  B  lie  on  the  same  streamline,   ^  =  +  f0   . 

A  '  B 

'xhe  only  way  in  which  these  two  conditions  on  f   can  be  satis- 
fied is  if  f,    =  -  V/-Q=  0   .   oince  the  stream  function  vanish- 
'A      B— 

es  at  the  st^gnati^n  points,  it  also  vanishes  at  all  other 
points  on  the  body,  t  ince  they  are  on  the  same  streamline. 

The  velocity  potential  <f   on  the  other  hand  varies  from 
a  maximum  negative  value  at  A  to  a  maximum  positive  value  at 
B.   i'he  value  fL   then  defines  the  amplitude  of  the  <f  func- 
tion on  the  contour  of  the  body,   i'his  amplitude  may  be 
arbitrarily  prescribed,  and  will  then  fix  the  size  of  the 
near  circle  contour  in  relation  to  the  fixed  spacing  J  -  2lf  . 

Let  K,  tXn.y-rj)  =  ii,D  and  so  forth.   Then,  at  the  point 
Lao  lx> 

B  from  the  foregoing  discussion; 


<PP  =  xK  cos/8  +  yR  sinp  +  AK1B  -  BJ1R  =  arbitrary  value 

A5a 


■B  =  x^   cosp  -i-  yB  bJ.xir    t  ^1B  -  «u 


y    =  -xB  sin^  +  yB  cos/3    +  AJlfi  +  BRlfi     =  0  A5b 

v     =  cos/3    +  Ah2B  -  BJpB     =  0  Ai?c 


v     =  sin/?  -  AJ2B  -  B^2B     =  0  A5d 
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i'he  last  two  of  the  above  equations  may  be  solved  simul- 
taneously for  the  unknown  constants  a  ana  B  with  the  results: 

R2B  cos/g  +  J2B  Sin^ 

A=  2 ~2 A6a 

*2B   +  J2B 

+  R2B  sin^3  +  j   cos£ 

B  =  ^ P A6b 

*2B     +  J2B 


substituting  these  results  into  tne  expressions  for  <#  and  V^ 
yields  two  equations  in  which  the  two  unknowns  are  the 
coordinates  j-q  ana  y  of  the  aft  stagnation  point.   The 
equations  arc  trancenaental,  however,  so  that  an  iteration 
technique  is  necessary  for  solution.   Care  was  necessary  in 
this  area  as  all  the  functions  are  multi-valued  in  y,  and  the 
result  could  converge  to  a  point  differing  by  some  integral 
multiple  of  2  Tf  in  the  y  dimension  from  the  intended  point, 
ihe  general  method  involves  starting  with  low  values  of 
in  effect  close  to  the  isolated  case  where  the  location  of 
the  stagnation  points  :s  known,  una  working  up  to  the  desired 
value  by  modest  increments. 

■Location  of  Point--  on  the  Contour 

^nce  the  location  of  the  aft  stagnation  point  has  been 
found,  the  constants  A  ana  B  can  be  evaluated  from  equations 
A6a  ana  A6b  .   i'he  other  points  on  the  contour  can  then  be 
located  from  equations  A5a  and  A5b  by  setting  ^=0  and 
assigning  a  suitable  value  to  ep  .   mis  computer  program 
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divides  the  interval  from  -f  to  +  \     into  20  equal  parts. 
ihis  yields  20  points  along  each  surface  of  the  airfoil,  or 
40  points  in  all.   As  pointed  out  in  the  text,  it  would  be 
helpful  to  subdivide  some  of  these  intervals  near  the  leading 
edge  to  get  a  better  check  on  the  maximum  velocity  in  this 
critical  region. 

I'he  equations  for  the  location  of  these  contuur  points 
are  again  trancendental,  and  a  similar  iterative  technique 
is  necessary  for  solution,   once  the  coordinates  of  a  point 
are  known,  the  velocity  components  are  found  from  equations 
A5c  and  A5d«   ihe  final  velocity  v  anu  the  slope  X  of  the 
contour  at  that  point  are  found  from  the  relations: 


v  =  f 


2     2 
v   +  v  A8a 

x     y 


v 
oC  =  arctan  -^  A8b 

x 

At  the  stagnation  points,  this  relation  becomes  indeterminate 

and  L'Hopital's  rule  is  used  for  a  solution. 

■Basic  relations  in  the  Circle  ilane 

it  is  well  known  that  the  flow  around  an  infinite  row 
of  arbitrary  but  identical  contours  can  in  principle  be 
mapped  into  a  corresponding  flow  about  a  single  unit  circle. 
if  the  flow  is  of  zero  circulation,  and  the  direction  of  flow 
far  from  the  contours  is  ^3  then  the  corresponding  flow  in 
the  circle  plane  is  expressed  b^  the  following  complex  po- 
tential and  velocity: 
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-  f  *  If-. ^ta^-J)  ♦..♦*"  inf^j 


« =  iv  v  ••*'  -  -^(j^  ♦  e+i%-^J  ^ 


Splitting  these  into  real  and  imaginary  parts,  and  placing  r 
=a  =  1  leads  to  the  result  that  along  the  unit  circle, 

Vu  =  0  AlOa 

AlOb 
f  =  constant 

This  verifies    that   the   circle   is   in  fact   a  streamline, 
further,   along  this   circle: 

f/3  n    /  cosh^-f  cos  &  \  /?.,,  /    sin  &  \         ,, 

=  cos/3  In/ 7-7: — i   +  smP    2  arc tan  /  — t-t-tt  J       Alia 
'         (    cosh/  -  cos<^  y  \  sinh^/ 


*->£ 


cosPcosh^   sini?  +  sin  /Q  sinh  <7    cos.  ^ 

coeh  2^  -  cos  2  ^  Alib 


Singularities  and  sjba^nation  points  in  the  Circle  Plane 

io  find  the  stagnation  point  Ht    set  V^  =  0   .   xhen 
equation  Allb  gives: 


tan  <£>  =  tan  /3  tanh  >*  a12 


Substituting  this  into  equation  AllA  produces  a  result  which 
can  be  reduced  to  the  form: 


62 


4B  =  coat    In  (i  !  °°^  )  +  einf.2  arctan(2fi?j  *l3a 


wnere 


I  =  ¥  £(  cosh  2^  +  cos  2/3)  A13b 

1'hese  equations  then  must  be  solved  for  T  then  f    for 
assigned  values  of  ft    and  <£L  which  must  be  the  same  as  the 
values  for  ft   and  <^L  previously  employed  in  the  near  circle 
plane.  Since  the  functions  are  trancendental ,  an  iterative 
technique  is  again  required.  Newton's  inetnod  is  employed, 
which  involves  calculation  of  tne  derivative  z — ■■  .  De- 
noting  this  by  the  symbol  QrR  then: 


B 


2 


B=   (  2?2  -  cos  2?)2   - 


A14 


Applying  the  equations  in  this  form  however,  requires  an 
excessive  number  of  significant  figures  in  the  evaluation  of 
T  .   io  circumvent  this  problem,  the  following  subltitution 
was  used  to  eliminate  1'. 


I  =  cos  fi   +^T  A15 


xhe  previous  equation  then  become 

r\  *    n     f ~2   cos  /3      +AT1        .     r>  0  .      /sin^ 

4>B  =  cos£    lnjf- jl -/+  Binf'2  arctanf- ^T+il  j| 


Al6a 
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<p'  .   /2  \/oosfl  +AT     )      (  cos  ■••  +  41-) >a6b 

xhe  problem  of  finding  I  for  assigned  values  of  /*  and<£ 
b„  Newton's  method  has  no  difficulty  with  regard  to  signifi- 
cant figures,   xhe  parameter/4  and  its  hyperbolic  functions 
may  now  be  found  from; 

sinh^  =  Yl'  (  2  cos/3  +^T)   =  ^  i   (  cosh  2/-1) 

A17 

xhe  initial  estimate  for  iiewton's  method  is  obtained  by 
neglecting  I  in  the  last  term  of  equation  Al6a  •  xhe 
resulting  equation  may  be  solved  explicitly  for  A  T  . 

^2  cos? 

(J  -1) 


where 


P  =  ($-2?%™?     . 


cos  ft 


■ 


ihese  expressions  provide  an  excellent  first  estimate  for 

T,  and  in  fact  for  the  special  case  otfis   0  ,  they  are 
exact. 

with  ft   and  ft*  known,  the  location  of  any  point  @. 
having  an  assigned  value  of  j.      is  defined  implicitly  by 
equation  Alia.   'Alo  establish  the  location  of  points  corres- 
ponding to  those  in  the  near  circle  plane,  all  that  is  neces- 
sary  is  to  have  the  respective  values  or  tf.    be  the  same  in 
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the  two  planes.  As  before,  an  iterative  solution  is  re- 
quired . 

Circulation 

'j;he  complex  potential  ana  velocity  for  a  purely  cir- 
culatory flow  of  unit  circulation  in  the  circle  plane  are 
as  folloxvs: 


a«  =  -  ^§sjl1h  / j ;  -  e^  j  il9a 


d^wj  -  2i  cos  fl    sinh  P  S 

d-5  54  -  2   cosh/*.  J2+l 


A19b 


i£> 

=  r  e  A19c 


Separating  these  into  real  and  imaginary  parts,  it  is  again 
found  that  the  unit  circle  r  =  1  is  a  streamline,  along 
which  the  velocity  distribution  is  given  by 

a  v      -  sinh  2#*  *2Q 

a  &~     cosh  2f    -   cos  2& 

This  circulatory  component  must  be  evaluated  at  the  same 
points  &.    previously  establishec^around  the  periphery  of  the 
circle,   iiy  superposition,  the  total  velocity  at  these  points 
is  then  given  by; 
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ota^nati.n  Points  with  circulation  in  uircle  ilane 

Just  as  the  stagnation  points  A  ana  b  are   defined  by 
the  condition  V   =  0  ,  so  are  the  stagnation  points  i,  and  T 
defined  b^  the  corresponding  condition  2£=  °  •   introducing 
this  condition  into  the  foregoing  equations  yields  the 
essential  result 

9B  ~  ^T  *  ^L  "  ^A  =  ^   =  arcsin/^°-  ■     Sinh  2^CQS4 

L-     cosh,^  cos  /3 

A22 


unce  the  points  L  and  !i  have  been  located  in  the  circle 
plane,  the  corresponding  values  of  <f  and  T;,.  m£*y  De  found 
from  equation  Alia,   introducing  these  \alues  into  the  near 
circle  plane  leads  to  the  location  of  the  points  L  and  1  in 
that  plane  by  solving  equations  Ma  and  A4b. 

Kaypin^  between  Circle  and  i\iear  Circle  Planes 

Happing'  between  the  circle  or  jf  plane  and  the  near 
circle  or  z  plane,  as  implied  in  the  preceeding  discussion, 
is  accomplished  despite  the  fact  that  no  explicit  function 
for  z  in  terms  of  K   or  vice-versa  is  known.   -»-he  mapping  is 
done  implicitly  through  the  fact  that  the  non-circulatory 
flows  in  the  two  planes  must  correspond,  or   in  other  words, 
through  the  condition  that 

vH  z)  =  w  (  $ )  A2pa 
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Different iating  this  gives 


idwi  =    Jml 

dzl  , dz  i 

to? 


A23b 


But 


and 


dwj  =  v  ^25c 

dz 


|M  1  =     V  A23d 

Jdj|  * 

Hence  the  mapping  derivative  becomes 

dz  |  =  Jjt  A2'3e 

d^  (    v 

xhe  flew  with  circulation  has  been  completely  defined 
in  the  foregoing  develoijment  for  the  circle  plane,   in  the 
near  circle  plane,  however,  only  the  non-circulatory  flow 
has  been  explicitly  defined,   as  exxjlained  in  the  text,  add- 
ing circulation  terms  to  the  complex  potential  in  the  near 
circle  plane  produces  unwanted  distortions.   Therefore,  the 
method  used  is  to  map  the  known  flow  with  circulation  in  the 
circle  plane  back  into  the  near  circle  plane  using  the  map- 
ping relations  previously  established.   These  mapping  rela- 
tions are  of  course  independent  of  circulation.   ±hus,  for 
the  final  flow  with  circulation 

v     U)  =  wn  (5)  =  w  (5)  +  B  Aw  (5)  A24a 
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Differentiating  this  gives 

—£  =  '^      o    d  a   ' 


dz  idz 

^3 


substituting  in  terms  of  velocities,  and  introducing  proper 
algebraic  signs  gives  finally  for  the  velocity  in  the  near 
circle  plane: 

1/   =  v*  ±  8  *v*       _  A„ /v  \ 

"    v#°      =  v  +  Bo^va(v^y        ^ 


Basic  delations  in  the  Cascade  Plane 

The  basic  idea  of  this  entire  method  is  to  map  the  near 

1 
circles  into  a  cascade  of  airfoils  in  much  the  same  manner 

that  the  Joukowski  transformation  may  be  used  to  map  a  single 

circle  into  an  airfoil,   ihe  equations  for  the  transformation 

are  identical  in  form  to  those  already  developed  for  the 

flow  in  the  near  circle  plane.  However,  all  coordinates  are 

evaluated  with  respect  to  axes  shifted  slightly  from  the 

original  position.  Denoting  these  shifted  coordinates  by 

primes,  the  mapping  function  employed  is  specifically: 


§*  =  f    +  aW   =   e"1^'   +  (A*  +  iB')  J^U1) 


A25a 


4^7  =  u  -  iu  =  e"1^'  +  (A1  +  iB')  F0(ss')         A25b 
dz'    x     y  2 
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where 


F1(zt)  =  1  coth  |-  A26a 

#2(.')  =  -*  each2  fl  A26fc 


Note  that  the  above  functions  are  exactly  of  the  same  form 
as  equations  a1  and  A2  .   All  subsequent  results  may  readily 
be  obtained  through  this  analogy. 

.Determination  of  Airfoil  ohape 

In  order  to  transform  the  near  circle  into  the  required 
airfoil  shape,  it  is  necessary  to  locate  the  singular  points 
P  and  \   of  the  mapping  transformation  in  the  proper  relation 
to  the  stagnation  points  Jj  and  i*  of  the  flow,   normally,  P 
is  locatea  a  short  distance  inside  the  contour,  just  opposite 
L  on  the  surface.   Similarly,  point  Q  lies  inside  the  con- 
tour on   a  line  normal  to  the  contour  at  i.  .aefer  to  i-igure 
5.   J-he  significance  of  t^ese  locations  may  be  qualitatively 
understood  from  the  fact  that  as  the  relative  distance  -— 
=»  £  is  made  small,  the  radius  of  curvature  at  the  nose  of 
the  airfoil  decreases  rapidly,  and  in  the  limit,  as  £  =  0 
tne  nose  becomes  cusped.   Similarly,  the  relative  distance 
r*  =  P.&       controls  trie  tail  curvature,   in  the  limit  as 
fi\  =   0   the  tail  becomes  cusped  in  the  manner  of  a  uoukow- 
ski  airfoil,   xhe  origin  of  coordinates  for  the  transforma- 
tion is  shifted  from  the  original  location,  the  midpoint 
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between  points  A  and  B  ,  to  the  new  location,  the  midpoint 
between  points  P  and  Q. 

Velocity  and  ^lope  ^long  the  Airfoil 

'ihe  complex  potential  including  circulation  ii-  the  cas- 
cade and  near  circle  planes  may  be  equated 

w   (t*)  =  w  (z)  A27a 

c  *  c 


Differentiating 

dw 

awc     =  ~df  A27b 

d$  *  /dS  »  j  /dz'  I 

Uz '     /Uz   / 

Le<t  (C¥c  \=  Ve"1  conjugate   of  velocity  A~7 

\dy*l  on  airfoil  ^'c        « 


\dz  /  conjugate  of  velocity 

on  near  circle 


-i *  A27e 

*.27f 


ii'rom  which,    alter  separating  amplitudes  ^nd  phase  angles: 


u 


o   =    cK^OC 


ii28b 
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APPMJDlA  a 

.Fortran  Program  "Cascade"  and  Operating  instructions 

on   the  following  pages  is  the  FQi&\tiAH  program  "Cascade" 

used  in  deriving  the  data  used  in  this  study,   xnput  data  is 
bv  means  of  data  cards  inserted  at  the  end  of  the  program. 
Any  number  of  data  cards  may  be  used,   -the  first  card  indica- 
tes how  many  additional  data  cards  are  to  be  read. 

iiaeh  additional  card  lists  the  number  of  values  of  the 
parameters   ,    ,  £     ,    and  iL     to  be  computed,  their 
initial  values,  and  the  increment  by  which  each  value  increa- 
ses.  Also  listed  are  the  print  control  quantities  NUPJfci  and 
NE&PR,  which  control  the  extent  of  the  print  out  as  indica- 
ted in  the  tabic  of  symbols,   ihe  compute*  then  calculates 
results  for  every  possible  combination  of  the  five  parameters, 
prints  the  results,  and  then  proceeds  to  the  next  card. 
ihe  standard  print  out  contains  the  five  inpit  parameters, 
the  solidity  and  stagger  angle  of  the  cascade,  the   airfoil 
thickness,  and  the  coordinates,  slope,  velocity,  and  boundary 
layer  parameter  for  each  of  the  42  points  of  ihe  airfoil. 

Condensed  Print  out 

when  surveying  large  numbers  of  cascades,  as  was  done 
for  this  study,  it  is  convenient  to  condense  the  data  to 
reduce  the  volume  of  output,   rage  82  shows  the  modifications 
made  to  the  basic  program  to  accomplish  this,   .briefly,  all 
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of  the  standard  print  commands  were  removed.   j.hen  commands 
to  print  the  heading  only  were  inserted  early  in  the  program 
right  after  changes  in  sigma.   The  modified  print  out  then 
contains  no  airfoil  data  at  all  except  for  the  maximum  value 
of  the  boundary  layer  parameter,   operation  of  the  program 
is  identical  with  the  basic  program. 

iterative  ocheme  for  Constant  xhickness  airfoils 

Page  84  shows  further  modifications  made  to  the  basic 
program  to  produce  cascades  with  the  desired  airfoil  thick- 
ness,  in  effect  this  eliminates  the  parameter  £  since  it 
is  varied  as  necessary  to  produce  the  required  thicknesses, 
l'he  initial  value  of  £  is  still  read  from  the  data  card,  but 
the  column  which  indicates  the  number  of  £,  to  be  calculated 
no   longer  has  significance,   l'he  program  as  written  will 
demand  thicknesses  from  4>  to  2A°/<>   in  4/"  increments,   uther 
values  for  thickness  can  be  obtained  by  modifying  the  appro- 
priate commands.  No  means  of  doing  this  by  data  cards  is 
provided.   The  print  oux  is  of  the  condensed  form  described 
previously . 
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.JOB  JENISTA 

PROGRAM  CASCAOE 

DIMENSION  XI(U2JfVIU2}.RI(U2),V3J«c2)tALI(U2),AHMU2).OELI(U2),TEI 

ligifin  u?rjc^iJuifI<4i5wE:£&5f?l5>f?h^lN^I,CU2,-v*Ifk2>  •vNnu 

201  READ  20U.IB 

204  FORMAT  (112) 

2021AKI?  g§^NB^|^§,NBZ'NE,>»NAK»BE1«SG1»Bi1tEPl,AKl.BEI,SGItBZI,EBI, 

205  FORMAT  tsf2,  1CF5. 2,212) 

PART  I A  CONSTANTS  FOR  NEAR  CIRCLE 
BEJ*BE1 

002  J-l.NBE 

BE*  (3.1U1593/180.)  •  BEJ 

SNBE*SINF(BE) 

CNBE*COSF(BE) 

SGK*SG1 

003  K*1,NSG 
PHEB*SGK#3,1U1593 
C22*SNBE 
C23*CNBE 

IF(K-I)  401, 401,402 
401  0RUP*1.5708 

DALP*0. 

SIP*0. 

RUP*0. 

ALP*BE 
U02    IF(SGK-SIP-.5)    U04,404,403 
U03   SIG*SIP*.5 

RU=RUP«-DRUP».5 

ALU*ALP*DALP«.5 

IF(NEXPR)  52C, 520,518 

518  WRITE  OUTPUT  TAPE  2,519, ORUP,OALP, SIP, SIG,RUP,RU, ALP .ALU 

519  FORMAT!/   7H  DRUP  *F8.4.   5X  6H0ALP  *  F9.4.  /   6H  SIP  «F  9.4,  5X 

1  5HSIG  *F10.4,  5X  5HRUP  *F10.4,  5X  4HRU  *F11.4,   5X  5HALP  -F10.4, 

2  5X  5HALU  *F10.4/   ) 

520  CALL  BAKER  ( SIG. XU,YU,RU ,ALU,AU,BU,R1U, AJ1U,R2U,AJ2U,. 000001, 
1    A2U,ITU,C22,C23,NEXPR) 

DRUP*(RU-RUP)/.5 

DALP*(ALU-ALP)/.5 

SIP*SIG 

RUP*RU 

ALP*ALU 

IF(NEXPR)  522,522,521 

521  WRITE  OUTPUT  TAPE  2, 519, ORUP,OALP,SIP, SIG, RUP.RU, ALP.ALU 

522  CONTINUE 
GO  TO  U02 

404  RB*RUP*DRUP«MSGK-SIP) 
ALB*ALP*DALPMSGK-SIP) 
IF(NEXPR)  525,525,523 

523  WRITE  OUTPUT  TAPE  2,52U, DRUP,DALP,SI P, SGK,RUP,RB, ALP, ALB 

524  FORMATS/   7H  DRUP  *F8.4,   5X  6HDALP  *  F9.4,  /   6H  SIP  «F  9.4,  5X 

1  5HSGK  *F10.4,  5X  5HRUP  *F10.4,  5X  4HRB  *F11.4,   5X  5HALP  *F10.4, 

2  5X  5HALB  *F10.4/   )  .     . 

525  CALL  BAKER  (     SGK,XB,YB,RB, ALB, A,B, RIB, AJ1 B,R2B,AJ2B, .O00CC1 , 
1    A2,IT,C22fC23,NEXPR) 

DRUP* ( RB-RUP ) / ( SGK-S I P ) 

DALP* ( ALB-ALP ) / ( SGK-SI P) 

SIP*SGK 

RUP*RB 

ALP*ALB 

IF(NEXPR)  527,527,526 

526  WRITE  OUTPUT  TAPE  2, 524, DRUP, OALP, SIP, SGK.RUP, RB, ALP ,ALB 

527  CONTINUE 
SNALB*YB/RB 
CNALB*XB/RB 

ALBO    »t180./3. 141593)*   ALB 

XA*-XB 

YA»-YB 
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ALAsALB*3.1U1593 

RA-RB 

XI(20)*XA 

XI(UO)=XB 

YI(2C)sYA 

YI(40)*YB 

RM20)»RA 

RI(UQ)«RB 

ALI(20)*ALA 

ALI (U0)=ALB^6. 283186 

VI(2C)»0. 

VI(UC)«C. 

R3B«-2.MR1B«R2B-AJ1B«AJ2B) 

AJ3B»-2,*(AJie«R2B*RlB»AJ2B) 

DVXB«A»R3B-B«AJ3B 

DVYB*B«R3B*A«AJ3B 

DVZB*SQRTF<DVXB»DVXB*DVYB»DVYB) 

AHI(l*0)*ASINF(DVYB/DVZB)-1#5707965 

AHI(20)«AHI(U0)^3.Uil593 

A22=A»A+B«B 

A23«RWRWAJ1«AJ1 

A2U  =  R2«R2<-AJ2»AJ2 

W2ZB*2.«SGRTF  ( A22*A23«A2U ) 

IF  (N0PR)120,120,121 

120  WRITE  OUTPUT  TAPE  2.  122,  BE J,SKBE,CNBE , SGK, PHE 6, A ,B ,XB,Y8,  Re, ALBC 
It  A2t  IT 

1220F0RMAT  (1H1///  46X  28HRESULTS  IN  NEAR  CIRCLE  PLANE  // 

1  10X  6HBETA  =F6.2,  5H  DEC  8X  10HSIN  BETA  *  F7.U,  8X 

2  10HC0S  BETA  *  F7.U,  8X  7HSIGMA  »  F7.3.  / 

3  10X  8HPHKB)  ■  F8.4,  9X  3HA  *  F1U.U.  8X  3HB  -  F1U.U,  // 
k       12X  1HI  9X  1HX  UX  1HY  llX  1HR  5X  12HLAMB0A,  DEC   6X  1 MV 

5  6X  11HALPHA.  DEC  3X5HERR0R  SX  10HITS   SHIFT  // 

6  12X  1H0  3F12.U,  F12.2,  6X  6H0.C000  8XUH»     E12.2,  16, 

7  6H      •) 

121  CONTINUE 
PART    IB    POINTS    ON   NEAR    CIRCLE 

DC    12    I«1.19 
IF(I-l)    13.13.1U 

13  PHN=.9«PHE8 
AL*ALB*.U5 
R»RB 
GC    TO    15 

14  AI»I 
PHN=PHEB»(1.-AI/10.) 
AL*AL-M  3. 141 593  +  ALB-AL)/ (21  .-AI ) 

15  X*R»COSF(AL) 

XALLSlsR(xiYfVX,VYfVtC22,C23,A,BtPHN,PHEBtITSfSNAL8fCNALBfUt 
1     ITSIDE,ERS,NEXPR) 
XKI  >=X 
YI(I)*Y 

R*SGRTF(X»X*V«Y) 
Rid  )«R 
VI(I)=V 

AHK  I)»ATANF(VY/VX) 
IF(VX)49,50.50       mt  mm^ 

49  AHK  I)*AHI(I)*3. 141593 

50  CONTINUE 
AHD  -M80./3. 141593)*         AHIC I) 
AL=ATANF(Y/X) 
IF(X)36.37,37 

36  AL*AL*3. 141593 

37  CONTINUE  „    ,.         ■_. 
ALD           *(180./3.1U1593)*  AL 
ALI(I)*AL 
XI(I*20)--(XI(I) ) 
YHI*20)«-(YI(I)) 
VIU*20)*-(VIII>),    ,,  __._ 
ALI(I*20)«ALI(I)*3.1U1593 
AHK  1*20)    *    AHK  I) 
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RI   (1*20)  *  RIII) 

IF  (NOPR)  12»4,12U,12  „ 

124  WRITE  OUTPUT  TAPE  2,  ^S.I.X.Y.R.ALH. V, AHD.ERS, ITS. I TSIOE 

125  FORMAT  (I13.2F12.Ut  F12.  2,  F12.M.F  J2  v2.E12.  2,  21  6) 
12  CONTINUE 

PART  1 1 A   CONSTANTS  FCR  CIRCLE 
IF(NEXPR)  530,530,528 

528  WRITE  OUTPUT  TAPE  2,  529 

529  FORMAT    (///) 

530  IT3»0 

POWER   *    <PHEe-2.«BE»SN8E 1/CNBE 
DET*(2.«CNBE)/(EXPF(P0WER)-1.) 

16  ARG    1    «    1.*(2.*CNBE)/0ET 

ARG2    »    SNBE/(CNBE^DET)  ma       mmhM 

PHE*    CNBE^LOGFt ARG  1  >*SNB£«2 . «AT ANF ( ARG2 > 

ER3*PHE-PHEB 

RER3*ABSF(ER3)/PHEB 

IT3MT3+1 

T«CNBE+OET 

Q1*CNBE*T 

Q2*l .♦Q1«DET 

PHP*(-2./DET)«(T/Q1)»(T/Q2) 

DETNU   *    DET-ER3    /PHP 

505  WRiTEXOUTPUT6TAPE,2?506,IT3,0ET.PHE.RER3,O|^^^      t¥Auoca,      -la 

506  F0RMAT(6H    IT3   *IU,5X5H0£T   -E10.U, 5X5HPHE   «F1 0.6.5X6HRER3   »E10.U, 
1       5X7HDETNU    *E1C.U) 

516    IFU000001-RER3)    17,19.19 

17  IFUT3-25)  18.19,19 

18  IF(OETNU)  60,61.61 

60  0ET*DET/2. 
GO  TO  16 

61  DETsOETNU 
GO  TO  16 

19  CH2G*1.+2.«Q1*0ET 

SH2G*SQRTF  (CH2G«CH2G-1 . ) 
SHG  *  SQRTF  ( (CH2G-1 . )/2. j 
CHG  »  SQRTF  ( (CH2G+1. )/2. ) 
C2=CH2G 
C3=SH2G 
CU'SHG 

TEB=ATANF((CU/C5)MC22/:23)1   tc 

TEBD    *<180./3.141593)«       TEB 

CNTEB=COSF(TEB) 

TEA=TEB*3.1W1593 

TEI(20)»TEA         „ft„,ax 

TEI (U0)*TEB*6. 283186 

VTB-O. 

DVB=-(C3>/(C2-CCSF(2.«TEB)) 

OVA=DVB 

VTI(20)*0. 

VTI(UO)*0. 

0VIC20)=0VA 

DVIUO)*DVB 

A22=C23»C5»CCSF|TEB 

A23=C22»CU»SINF(TEB) 

A24=C2-COSF(2.«TEB> 

W2TB=U.»(A22*A23)/A2U 

1260WRITE°0UTPUT6fAP^2ri27,SHG,CHG,DET,PHP.TEBp,DVB.^ 
270FCRMAT(//UlX37HCORRESPON6lNG RESULTS    IN   tlRCL|    PLANE// 

1       10X    12HSINH    GAMMA   *    F8.U,    5X    12HC0SH    GAMMA    «    F8.U,    5X    7HDEL    T 

5      E19.2,     113) 
128    CONTINUE         _,   „M    .,__,  c 
C      PART    I  IB      POINTS   ON   CIRCLE 


75 


DO    20    1*1,19 

AIM 

PHN  =  PHEBM 1.-AI/10.) 

IF(I-1)21,21,22 

21  TEOLD   =    TEB 

TE    *    TEB    ♦    SCRTF    ( .2»PHEB/W2TB) 
GO    TO   23 

22  TEOLD   «    TE 

TENU    1    =    TE-.1»PHEB/VTE 
M    IF    (TEA-TENU    1)    53,53,51* 

53  TE    =    (TE*TEA)/2. 
GO   TO   23 

54  TE    *    TENU    1 

23  ITU*0 

25  C6*SINF(TE) 
C7*C0SF(TE) 
C8«C0SF(2.«TE) 
A2»L0GF((C5*C7)/(C5-C7n 
A3*2.»ATANF(C6/C4) 
PH«C23»A2+C22»A3 
A2=-(C23«C5»C6)*C22»CU»:7 
VTE=U.»A2/(C2-C8) 
DPH=PH-PHN 

ITUMTU+1 
A2*ABSF(DPH)/PHEB 
IF(NEXPR)    511,511,512 

512  WRITE   OUTPUT    TAPE    2,    513,    I,    ITU,    TE»DPH,VTE.A2 

513  FORMAT    (       6H         I    *I4,    5X5HITU    «l4,    6XWHTE   »E10. 
1       5X5HVTE    *£10.U,    5XUHA2    *E10.4) 

511  IFUTU-25)  2U, 111.111 

24  IF  (  ,000001  -  A2)  26,  Ml,  111 

26  TENU  *  TE-DPK/VTE 

IF    (TENU  -    TEOLD    )    51,51  ,55 
51    TE    *    (TE*TE0L0)/2. 
GO    TO   25 

55  IF    (TEA-TENU)    56,56,57 

56  TE    *    (TE*TEA)/2. 
GO    TO   25 

57  TE    •    TENU 
GO   TO   25 

HI    TED  =(180./3.141593)»      TE 

TEI(I)    =    TE 
VTK  I)*VTE 
DVK  I)=-C3/(C2-C8) 
TEI(I*20)*TEI(I)*3.1U1593 
VTK  K20)*-(VTKI)) 
DVK  K20)    =    CVKI) 
IF    (NOPR)    129,129,20 

129  WRITE   OUTPUT    TAPE    2,    130,  Iv  TEDf  VTE.DVK  D.A2.ITU 

130  FORMAT    (116,    F17.2,    F21.4,    F20.4,    E19.2,    113) 
20    CONTINUE 

PARTIIIA      CONSTANTS    FOR    AIRFOIL 

BZL=BZ1 

D029L=1,NBZ 

BAM=ABSF(BZL) 

IF(8AM-. 001)65, 85, 86 
85    TAU=0. 

XT*XB 

YT*YB 

XI(42)=XT 

YI(42)*YT 

VT*0. 

AHT*AHI(40) 

AHL-AHI (20) 

AHI(41)*AHL 

AHK42)*AHT 

TET=TEB*6. 283186 

PHT*PHEB 

XL*XA 

YL»YA 
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U,    6X5HDPH    »E10.4, 


XICU1 )»XL 

Yl(41)=YL 

VL  =  0. 
TEL=TEA 
PHL=-PHEB 
GO  TO  87 
86  CONTINUE 
C2=CH2G 
C3=SH2G 
CUsSHG 
C5=ChG 
C22=SN8E 
C23-CNBE 

TAy  =  ASINF((BZL/U.M(SH2G/CHG)«CCNTEB/CNBe)l 

TET=TEB-TAU*6«28  3186 

TEL=TEA*TAU 

A2=SINF(TET) 

A3*COSF(TET) 

Ai*=(C5*A3)/(C5-A3) 

A5-A2/CU 

PHT=C23«LOGF(A4)*C22»2.»ATANF(A5) 

A2=SINFCTEL) 

A3=COSF(TEL) 

AU=(C5*A3)/(C5-A3) 

A5=A2/C4 

PHL=C23»LOGF(AU)+C22»2.«ATANF(A5) 

ENT=30.*10.»PHT/PHEB 

1  =  39 

32  AI  =  I 

IF(ENT-AI)  3C.31.31 

30  1*1-1 

GO  TO  32 

31  A2*(TET-TEI(I))/(TEI  (1  +  1  l-TEK  X )  ) 
RT-RI(I)+A2«(RI(I+1)-RI(  I)} 
ALT*ALim+A2«(ALX(I*1)-AU<m 
XT=RT«COSF(ALT) 
YT=RT«SINFCALT) 
ENL=30.*10.»PHL/PHEB 

1=21 

33  AI«I 
IF(AI-ENL)3U,35.35 

34  1=1+1 
G0T033 

35  A2=(TEL-TEI(  X))/(TEI(I+1  1— TEI  fill 
RL»RI(I >*A2»(RI< I*1)-RI(  I)) 
ALL=ALI (I>*A2«(ALI(I*l)-ALI(in 
XL=RL»COSF(ALL) 
YL=RL»SINF(ALL) 

CALL  SSR<XT,YTtVX,VYtVT.C22.C23tA,B,PHT.PHEBtITStSNAlB.CNAL8.-1., 
1   ITSIDE,ERST,NEXPR) 
XI(42)=XT 
YI(    U2)=  YT 
AHT=ATANF(VY/VX) 
IF(VX)63. 64,64 

63  AHT=AHT*3. 141593 

64  AHI(U2)=AHT 
CALLSSR(XL.YL.VX,VY.VLfC22.C23.A,BfPHL.PHEB. ITS fSNAlB,CNALB,-1 .» 

1   ITSIOE.ERSL.NEXPR) 
XI(U1)=XL 
YI(41)*YL 
AHL=ATANF(VY/VX) 


IF(VX)65,66,66 
65  AHL=AHL*3. 141593 


66  AHI(U1)=AHL 
87  CONTINUE 

TNBE=SNBE/CNBE 

TNBE2-TNBE-B2L/2.  mm    m   mmm% 

CV2=( 1.*TNBE*TNBE)/( 1.*T NBE2»TNBE2) 

CV=SQRTF(CV2) 

EPM*EP1 
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DO    40    M*1,NEP 

AKN^AKl 

DC    41    N=1,NAK 

IF(EPM)    99,96,99 

98  IF(N-I)    99,99,41 

99  CONTINUE 
XP=XL-EPM»RB«SINF(AHL) 
YP*YL+EPM«RB«COSF(AHU 
XC»XT-AKN»EPK*RB«SINF(AHT) 
YQ=YT+AKN»EPM«RB«COSF(AHT) 
XZ=(XP*XQ)/2. 
YZ=<YP*YQ)/2. 
XGP=XQ-X2 

YQP*YC-YZ 

X*XQP 

Y  =  YCP 

C2=SINH(X) 

C3=C0SH(X) 

C6=SINF(Y) 

C7=C0SF(Y) 

Rl=(C2/(C3-C7))/2. 

AJ1={C6/(C3-C7))/(-2«) 

PX=C3»C7-1. 

PY=C2*C6 

PZ=PX«PX*PY»PY 

R2=PX/<PZM-2,)) 

AJ2=PY/(PZ«2.) 

QZ=R2*R2+AJ2«AJ2 

PU»Y*QZ-R2»AJUAJ2«R1 

PV=X»QZ-AJ2»AJl-R2«R1 

BEP-ATANF(PU/PV) 

C22=SINF(BEP) 

C23=COSF(BEP) 

QX=AJ2«C22-R2«C23 

QY=R2»C22*AJ2»C23 

APsQX/QZ 

BP=QY/QZ 

AP22=AP«AP4-BP»BP 

AP23=R1»R1*AJ1»AJ1 

AP24=R2«R2*AJ2«AJ2 

WP2Z8=2.»SCRTF(AP22«AP23«AP2U) 
PARTIIIB   POINTS  ON  AIRFOIL 

1  =  1 
42    X=XI (I)-XZ 

Y=YI(I)-YZ 

C2=SINH(X) 

C3=CCSH(X) 

C6=SINF(Y) 

C7=C0SF(Y) 

R1=C2/(2.»(C3-C7)) 

AJl=C6/(-2.»(C3-C7)l 

PX=C3*C7-1. 

PY=C2»C6 

PZ=PX»PX*PY«PY 

R2=PX/(PZ«<-2.)) 

AJ2=PY/(2.»PZ) 

XI I  (  I)=X»C23-»Y«C22  +  AP»R1-BP«AJ1 

ETI  (  I  )=-X»C22+Y*C23*AP»AJWBP»Rl 

A2=C23+AP»R2-BP*AJ2 

A3=C22-AP«AJ2-BP*R2 

A4=SQRTF( A2»42+A3«A3) 

IF(AU-.OOOI)     171,171,172 
171    RP3=-2.«(R1*R2-AJ1»AJ2) 

AJP3=-2.«(  AJl»R2«-Rl»AJ2) 

DUXP*AP«RP3-eP«AJP3 

DUYP=BP»RP3**P*AJP3 

DUP=SGRTF(DUXP»CUXP*DUYP«DUYP)  „ 

EPI ( I)=ASINF(DUYP/OUP)-1 . 5707965 

IF(ALI(n-ALI(30))  173,173tH3 
173  EPI(I)*EPI(l>*3. 141593 
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172 

600 
113 


88 

89 

38 
39 

U3 


59 
62 

U6 

47 


710 

716 
715 


GO    TO    113 

EPK  I)=ASINF    (  A3/AU) 
IF(A2)    600,113,113 
EPI(I)  =  3.H4l593-EPI(I) 
DELI(I)=(180./3.1U1593)» 
A6=ABSF(VTI( I)) 
IF(A6-.0001)    88,88,89 
A5»SQRTF(W2Ze/k2TB) 
GO    TO    38 
A5«VI(I )/VTI(I) 
IF(ABSFIAU)-.OOOI)    39,39 
IFtBAH-.OOOl )    43, 43, 44 

yyi(I>-SQRTF{W2ZB/WP2ZB) 
GO    TO    62 

R3=-2.»<R1«R2-AJ1»AJ2) 

AJ3=-2.»(AJ1 «R2*R1«AJ2) 

WP2Z=SQRTF<AP22»(R3«R3*A 

SHX=SINH(XI{ I )) 

CHX=CCSH(XI( I)) 

SNY=SINF(YI( I) ) 

CNY=COSF(YI( I)) 

Rl«(SHX/2.)/(CHX-CNY) 

AJl=(-SNY/2. )/(CHX-CNY) 

PX=CHX«CNY-1. 

PY=SHX«SNY 

PZ*PX»PX+PY«PY 

R2=PX/(-2.»PZ) 

AJ2=PY/(2.«PZ) 

R3*-2.»(R1«R2-AJ1»AJ2) 

AJ3=-2.»{ AJ1«R2*R1«AJ2) 

W2Z=SQRTF( (A*A*B«B)«(R3» 

SNTE=SINF(TEI(I) ) 

CNTE=COSF(TEI(I) ) 

SN2TE  =  SINF(2.»TEIU)) 

CN2TE=C0SF(2.»TEI(I)) 

B22=CNBE«CHG«CNTE 

B23=SNBE*SHG*SNTE 

B24=CH2G-CN2TE 

W2T=-4.«(B22+B23)/B24 

DW2T=(2.«SH2G»SN2TE)/(B2 

A25=W2Z-A5«A5»W2T 

A26=(DVIU)*A25)/A6-A5«A 

VVK  I)=(W2Z+8ZL«A26)/WP2 

GO   TO   62 

VVI ( I )  =  ( VI ( I  )+BZL»DVHI) 

IF(  1-42)    U6,»47,i*7 

1  =  1*1 

G0TC42 

CONTINUE 

DELX=XII(42)-XII (41) 

DELY=ETI(U2)-ETI(l4l) 

CHORD=SQRTF(CELX«OELX*DE 

SOLID=CHORD/6.283186 

SNA=OELY/CHORD 

CNA=CELX/CHOPD 

AROT=( 180./3.141593)»ASI 

STAG=M180./3.141593)»(BE 

AHTD  =  ( 18C./3.1U1S93) 

TETD  =  (  18C./3. 141593) 

AHLO  *(18C./3. 141593) 

TELO  =(18C./3.1M1593) 

TAUD  =  (  180./3.1M1593) 

BEPD  =  (  180. /3. 141593) 

IAA=    20 

00    717    IDD»1,    1C 

IAA    =IAA-1 

Is    35 

TOD-TANF  <DELI(IAA)«3.14 
AA*  XNHI  )-XNI(IAA)  -  (Y 
IF    (AA)    714,71*1,713 


(AHim-EPicm 


59 


J3«AJ3)) 


R3*AJ3»AJ3)) 


4»B24) 

5«CW2T 
Z 

•A5)/A4 


LY«DELY) 


NF    (SNA) 
P+ASINF(SNA)) 

•  AHT 

•  TET 

•  AHL 

•  TEL 

•  TAU 

•  BEP 


1593/180.) 

NIC  IAA)-YNICI))«TOD 
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jl .. 


■ 


■ 


713    1  =  1-1 

GC    TC   715 
71U    BBE=ATANF    (IYNKI  +  1)    -YN I ( I ) ) /(XNI ( I *1 )-XNI< I ) ) ) 

DOEL    *    .5«(BeE*OELI(IAA)»3. 141593/180.) 
l    ^    CC    =    ABSF((YNI(IAA)-YNI(  1*1 ) )/COSF t OCEL) ) 

717  DDKICO)    =    CC*    (XNKI  +  D-    XNI    (IAA)    -CYNIUAA)    -    YN I  (  I  ♦  1 )  )  • 
1  TANF(ODEL))»SINF    (BBE     ) 

WAXIE  =  OOI ( 1) 

DO  718    100  =  2,10 

718  WAXIE  ~    MAX1F  (WAXIE,  001(100)) 

"    WRITE  OUTPUT  TAPE  2,302,  BEJ,SGK,B2LfEPM, AKN , STAG, SOL IO,W AXI E 
302  FORMAT  ( 1  Hi //52X 15HC ASCAOE  RESULTS// 

1  18H  DESIGN  PARAMETERS/ 

2  9H    BETA  -F6.2,  5H  DEG.  6X7HSIGMA  «  F6.3,10X8M8  ZERO  ■ 

3  F6.2,  9X  9HEPSIL0N  *  F5.3,9X7HKAPPA  «  F6.3,// 
U    17H  OTHER  PARAMETERS/ 

5     11H    STAGGER«F6.2,5H  OEG.  4X  10HS0LI0ITY  -F6.3  ,7X17HMAX.  THI 

6CKNESS  =  F5.3) 
IF  (NOPR)  lliCtlUOflUI 
140  WRITE  OUTPUT  TAPE  2, 304, XB.YB.PHEB, XT, YT. AHTO, TET0,PHT,XL ,YLt AHLO, 

1   TELD.PHL^P.YP.XQ.YQ.TAUO^XZ.YZ.AP.BP.BEPD 
3040F0RMAT  (3XUHXB  =  F8.4  ,  1 1X«4HYB  *F8.4,  11X8HPHKB)  *F8.4./3X 

1  UHXT  =  F8.4.11X4HYT  *F8.U,  1 1 XI OHALPHAt T)  «  F6.2.  5H  DEG. 

2  2X10HTHETA(T)  =  F6.2,  5H  DEG.  2X8HPHKT)  *  F8.4./3X 

3  4HXL  =  F8.U.  11XUHYL  »  F8.4,  1 1X10HALPHAI L)  «  F6.2,  5H  OEG. 
U    2X10HTHETML)  =  F6.2,  5H  OEG.  2X  8HPHKL)  *  F8.H,/3X 

5  4HXP  =F8.4.11X4HYP  =F 8.4 , 1 1X4HXQ  *F8.4, 1 1X4HYQ  «F8.4,11X 

6  5HTAU  =F13.2.  5H  0EG.,/3X 

7  4HXZ  =  F8.4,  11XUHYZ  *  F8.4,  UXUHAP  *  F8.U,  11X4HBP  «  F8.U.11X 

8  12HBETA  PRIME  =  F6. 2,  5H  DEG.) 
1U1  WRITE  OUTPUT  TAPE  2,303 

3030F0RMAT  (//M2X  13HAIRF0IL  SHAPE  32X21HVEL0CITY  DISTRIBUTION/ 

1  88X  19HAT  DESIGN  CONDITION/  m M    . 

2  10X1HI17X1HX19X1HYUX1  1HSLCPE,  DEG.  11X6HC V/V2)  12X9HC V/ V2)»«2//) 
1=1 

DO  90  1=1,42 

PRX=(XII( I )-XII(41 ))/CHORO 
PRY=(ETI<  I)-ETl(UI))  /CHORD 
XNK  I)=PRX*CNA+PRY»SNA 
YNI (  I  )=PRY*CNA-PRX«SNA 
DELI(I)=DELI (I)-AROT 
VVI (I )=CV»VVI( I ) 
V2I (  I)  =  VVI(I)»VVI(I) 
90  CONTINUE 

WRITE0UTPUTTAPE2,305.(I.XNIC IltYNIC I)90ELICI)tVVICI}fV2I( lit 

305  FORMAT  (111,  2F20.4.  F20.2t  2F20.4) 
PARTIV   PARAMETER  CHANGES 
48  AKN=AKN+AKI 
44  CONTINUE 

EPM=EPM+EBI 
UO  CONTINUE 

BZL=BZl+BZI 
29  CONTINUE 

SGK=SGK+SGI 
3  CONTINUE 

BEJ=BEJ+BEI 
2  CONTINUE 

10=10*1  m    m 

IF(IB-ID)  203,203,202 

2PARTNV   ITERATION  SUBROUTINE  FOR  POINTS  ON  NEAR  CIRCLE 

SUBROUTINE  SSR  ( X, Y, VX,V Y, V.C22 ,C23, A, B, PHN, PHEB, ITS, SNALB,CNALB, 
1SI0E,ITSIDE,A2,  NEXPR) 
ITS  =  C 
ITSIDE=0 
IF(NEXPR)  8,8,517 

US  F0RMAT?//5HT  ITsMlHX  9X1HY  9X2HR1  7X3HAJ1  8X2HR2  7X3HAJ28X1HA 
1   9X1HB  8X2HVX   8X2HVY   8X2HA2   5X4HSIDE/) 
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- 


8  C2=SINH(X) 
C3=C0SH{X) 
C6=SINF(Y) 
C7=C0SF(Y> 
A2=C3-C7 
R1*C2/(A2»2.  ) 
AJl=C6/(A2*{-2.>) 
PX*C3»C7-1. 
PY=C2«C6 
PZ=PX*PX+PY*PY 
R2*PX/(-2.»PZ> 
AJ2=PY/t2.»PZ) 
PH=X«C23+Y»C22*A«Rl-B«AJl 
PS»-X»C22  +  Y«C23*A«AJHB»Rl 
DPH=PH-PHN 
VX=C23*A«R2-E*AJ2 
VY*C22-A«AJ2-B«R2 
V=-SIOE«SQRTF(VX»VX*-VY«VY) 
A2=( ABSF(DPH)*ABSF<PS))/PHEB 
ITS=ITS-H 

IF(NEXPR)    507,507,508 

l§irc^lT0Viri5FA^.2:5?!5!IJ'?5T'R1'*J,-R2'AJ2-*'e'wx'v¥'A2',TS,DE 

507    IFUTS-25)    9,100,100 

9  IF    (.000001-A2)     10,100.100 

10  DX=( VX*0PH-VY«PS)/(V»V) 
DYM VY*DPH+VX»PS)/(V«V) 
X=X-DX 

Y=Y-0Y 
GO    TO    8 

100  YY=Y«CNALB-X«SNALB 

IF    (YY*SIDE) 101,11,11 

101  ITSIDE=ITSIDE+1 
IFUTSIOE-3)     102,11,11 

102  YY=-YY 
XX=X»CNALB+Y*SNALB 
X=XX«CNALB-YY*SNAIB 
Y=YY»CNALB+XX«SNALB 
GC  TC  8 

11  RETURN 
ENO 

PART    VI     ITERATION    SUBROUTINE    FOR    INITIAL    STAGNATION    POINTS   ON    NEAR    CIRCH 
SUBROUTINE    BAKER     (  SGK  ,XB,YB,RB, ALB,A,B,R 1 , AJ 1 ,R2, AJ2, ERM , A2, 

1  IT,C22,C23,NEXPR) 

IF(NEXPR)515,515,500 

500  WRITE    OUTPUT    TAPE    2,501 

501  F0RMAT(//5H         IT    UX2HXB    8X2HYB    8X2HR1    7X3HAJ 1    8X2HR2    7X3HAJ2    8X1HA 
1       9X1HB    9X2HPS    8X2H0P    8X2HA2/) 

515    CONTINUE  , 
IT=0 

PHEB=SGK*3.U1593 
SNALB=    SINF(ALB) 
CNALB=    COSF(ALB) 
XB=RB*CNALB 
YB=RB»SNALB 
4    C2=SINH(XB) 
C3=C0SH(XB) 
C6=SINF(YB) 
C7=C0SF(YB) 
A2=C3-C7 
R1=C2    /(A2»2.) 
AJl=C6/(A2»(-2.) ) 
PX=C3«C7-1. 
PY=C2*C6 
PZ=PX*PX+PY«PY 
R2=PX/(PZ*<-2.)) 
AJ2=PY/{PZ*2.) 
QX*AJ2»C22-R2«C23 
QY*AJ2»C23*R2»C22 
QZ*R2*R2+AJ2»AJ2 
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A=QX/OZ 

B=QY/QZ 

PH=XB»C23+YB*C22+A«R1-B»AJ1 

PS=-XB»C22*YB»C23*A»AJHB«R1 

DP=PH-PHEB 

A2=( ABSF(DP)+ABSF(PS))/PHEB 

IT=IT+1 

IF(NEXPR)    5H4,51U,503 
503    WRITE    OUTPUT    TAPE    2, SOU, IT, XB, YB,R1 • AJ1 , R2, AJ2t A,B»PS,OPt A2 
SOU    FORMAT    (I5.1CF10.UtlE10.il) 
51U    CONTINUE 

IFUT-25)    5,7,7 

5  IF(ERM-A2)    6,7,7 

6  A2*C3-C7 
R1X=(1.-C3»C7)/(A2«A2»2.  ) 
RlY=(C2»C6)/(A2«A2M-2.) ) 
AJ1X=-R1Y 

AJ1Y=R1X 
PXX=C2«C7 
PXY=-C3«C6 
PYX=-PXY 
PYY=PXX 

PZX=2.«(PX»PXX+PY»PYX) 
PZY=2.»(PX«PXY*PY«PYY) 
R2X=(PZ»PXX-PX»PZX)/((PZ«PZ)«<-2.)) 
R2Y=(PZ«PXY-PX»PZY)/(PZ«PZ»(-2. )) 
AJ2X=-R2Y 
AJ2Y=R2X 

QXX=-R2X»C23+AJ2X»C22 
QXY=-R2Y»C23*AJ2Y»C22 
QYX=QXY 
QYY=-QXX 

0ZX=2.*(R2»R2X+AJ2»AJ2X) 
QZY=2.*(R2«R2Y+AJ2»AJ2Y) 
AX=(QZ»QXX-QX«QZX)/(GZ»QZ) 
AY=(QZ«QXY-QX*QZY)/(CZ«QZ) 
BX={QZ»QYX-OY»QZX)/(QZ»QZ) 
BY=(QZ*QYY-QY«CZY)/(QZ*QZ) 
PHX=C23+A»R1X*AX«R1-B«AJ 1X-8X«AJ1 
PHY=C22+A»R1Y*AY«R1-B»AJ1Y-BY»AJ1 
PSX=-C22+A*AJ1X+AX»AJ1+B»R1X*BX«R1 
PSY=C23*A»AJ  mAY»AJl*B»RlY*BY«Rl 
A2=PHX*PSY-PSX«PHY 
DX=(PSY«DP  -PHY«PS)/A2 
DY*(PHX»PS-PSX«DP  )/A2 
XB=XB-DX 
YB=YB-DY 
VG0  TO  U 

7  RE=SQRTF(XB»XB+YB«YB) 
ALB=ATANF< YB/XB) 
IF(XB)28,U5,M5 

28  ALB=ALB+3.1U1593 
U5  CONTINUE 

RETURN 

END 

END 


1  5  6  1  5  U5.  2.5 


5   .001   .01 


•  U 


15   0    .245   1  0 


..END 
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MODIFICATIONS   i-'oii  DOJiD.fciMS.bUJ  i^uiA*!   uUT 


PART  I A  CONSTANTS  FOR  NEAR  CIRCLE 
BEJ-BE1 

002  J=1,NBE 

BE*  (3.1U1593/180.)  •  BEJ 

SNBE=SINF(BEI 

CNBE=COSF<BE) 

SGK=SG1 

003  K=1,NSG 

740  WRITE  OUTPUT  TAPE  2.  741,  BEJ.  SGK 

741  FORMAT  (  1H1  /26HCASCADE  RESULTS  FOR  BETA  =»F6.2,5H  DEG.UX7HS IGMA  « 
1F6.39//75H    B  ZERO      KAPPA      EPSILON      SOLIDITY      THICKNE 
2SS      BL  PARAMETER/ 1 


DO  718    IDD  *  2f 10 
718  WAXIE  *  MAX1F  (WAXIEt  DDKIDD)) 


751  PRAM»  V2I(1) 

760  PRAM  *  MAX1F  (PRAM.  V2IJI))  n   u.wIF-  PRAM 

WRITE  OUTPUT  TAPE  2t    I55t,BZLt  AKNt  EPM,  SOLIDt  WAXIEt  PRAM 
755  FORMAT  ( 3X.F6.2.2F 10. 3t2F14.3tF16.%) 
PARTIV   PARAMETER  CHANGES 


HuJDlJj'IOATiuJtfS   tfvLi   Uuj.«.jiia.,x    iida-Ciii^iiSS 


87   CONTINUE 

TNBE*SNBE/CNBE 
TNBE2*TNBE-BZL/2. 

CV2*( 1.4TNBE»TNBE )/( 1*+TN6E2*TNBE2> 
CV*SQRTF(CV2) 
.    AKN=AK1 

DO    U1    N=1,NAK 

WAXIT    =    0 

THIK    =    .04 

EPM»EP1 

IF(EPM)    99,98,99 

98  IF(N-I)    99,99,41 

99  CONTINUE 


1  TANF(DDEL) )»SINF    ( BBE    ) 

WAXIE    =    DDK1) 
mmm    DO    718         IDO    =   2,  10 
718    WAXIE    =    MAX1F    (WAXIE,    D0ICID0I1 

756  IF  (WAXSG)  757,  757.  751 

757  IF  (THIK  -  WAXIE)  761,  751,  763 
761  IF  (WAXIT)  764,  764,  765 

764  EPM  =  EPM  -  .01 
GO  TO  99 

765  SfScS  EPM,^  -01»(  (WAXIE  -  THIK)/(WAXI£  -  WAXIT)) 
WAXSG  *  100.0 

GO  TO  99 
763  EPM  *  EPM  ♦  .01 

WAXIT  -  WAXIE 

GO  TO  99 
751  PRAM=  V2I (1) 

DO  760  1=  1,  19 
760  PRAM  =  MAX1F  (PRAM.  V2KI)) 

7.;.;  SSi^T°VIEULTJPicfA  I5§£,BZk'cAfNf4EPM'  S0LI0»  WAXIE,  PRAM 
755  FORMAT  ( 3X,F6.2,2F10. 3,2F14»3,F16.4) 

EPMsEPM+EBI 

THIK  -    THIK  ♦  .04 

WAXIT  '    0 

WAXSG  =  0 

IF  (THIK  -  .25)  99,  99,  48 

PARTIV   PARAMETER  CHANGES 

48  AKN=AKN+AKI 

41  CONTINUE 

BZL»BZL+BZI 

29  CONTINUE 

SGK=SGK+SGI 

3  CONTINUE 

BEJ=BEJ+BEI 

2  CONTINUE 

ID*ID*1 

„  ,  IF( IB-ID)  203,203,202 

203.  END 
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